Spin electric effects in molecular antiferromagnets 
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Molecular nanomagnets show clear signatures of coherent behavior and have a wide variety of 
effective low-energy spin Hamiltonians suitable for encoding qubits and implementing spin-based 
quantum information processing. At the nanoscale, the preferred mechanism for control of quantum 
systems is through application of electric fields, which are strong, can be locally applied, and rapidly 
switched. In this work, we provide the theoretical tools for the search for single molecule magnets 
suitable for electric control. By group-theoretical symmetry analysis we find that the spin-electric 
coupling in triangular molecules is governed by the modification of the exchange interaction, and 
is possible even in the absence of spin-orbit coupling. In pentagonal molecules the spin-electric 
coupling can exist only in the presence of spin-orbit interaction. This kind of coupling is allowed 
for both s — 1/2 and s = 3/2 spins at the magnetic centers. Within the Hubbard model, we 
find a relation between the spin-electric coupling and the properties of the chemical bonds in a 
molecule, suggesting that the best candidates for strong spin-electric coupling are molecules with 
nearly degenerate bond orbitals. We also investigate the possible experimental signatures of spin- 
electric coupling in nuclear magnetic resonance and electron spin resonance spectroscopy, as well 
as in the thermodynamic measurements of magnetization, electric polarization, and specific heat of 
the molecules. 

PACS numbers: 75.50.Xx, 03.67.Lx 



I. INTRODUCTION 

The control of coherent quantum dynamics is a nec- 
essary prerequisite for quantum information processing. 
This kind of control is achieved through coupling of the 
internal quantum degrees of freedom of a suitable micro- 
or mesoscopic system to an external classical or quantum 
field that can readily be manipulated on the characteris- 
tic spatial and temporal scales of the quantum system. 

The molecular nanomag nets (MNs) 0, Q represent a 
class of systems that show rich quantum behavior. At low 
energies, the MNs behave as a large spin or a system of 
only few interacting spins. The behavior of this spin sys- 
tem can be designed to some degree by altering the chem- 
ical structure of the molecules, and ranges from a single 
large spin with high anisotropy barrier, to small collec- 
tions of ferro- or antiferromagnetically coupled spins with 
various geometries and magnetic anisotropies. This ver- 
satility of available effective spin systems makes the MNs 
promising carriers of quantum information [3] . While the 
interaction with magnetic fields provides a straightfor- 
ward access to the spins in an MN, it is preferable to use 
electric fields for the quantum control of spins, since the 
electric fields are easier to control on the required short 
spatial and temporal scales. In this work, we explore the 
mechanisms of spin-electric coupling and study the ways 
in which an MN with strong spin-electric coupling can be 
identified. 

Quantum behavior of MNs is clearly manifested in the 
quantum tunneling of magnetization j44lll| . A prototyp- 
ical example of quantum tunneling of magnetization is 
the hysteresis loop of an MN with a large spin and high 
anisotropy barrier. The height of the barrier separating 



the degenerate states of different magnetization leads to 
long-lived spin configurations with nonzero magnetic mo- 
ment in the absence of external fields. The transitions 
between magnetization states in the MN driven through 
a hysteresis loop occur in tunneling events that involve 
coherent change of a many-spin state. These transitions 
have been observed as step-wise changes in magnetiza- 
tion in single-molecule ferromagnets 0, H Il2l - fl4j . Sim- 
ilar tunneling between spin configurations are predicted 
in antiferromagnetic molecules |15l |16| , and the observed 
hysteresis was explained in terms of the photon bottle- 
neck and Landau-Zener transitions fl7H2C| . The transi- 
tions between spin states are coherent processes and show 
the sig natures of interference between transition paths 
|2lj - j23| , as well as the effects of Berry phase in tunneling 

Spin systems within molecular nanomagnets offer a 
number of attractive features for studying the quantum 
coherence and for the applications in quantum informa- 
tion processing [3j . A wide variety of spin states and cou- 
plings between them allows for encoding qubits. Chem- 
ical manipulation offers a way to modify the structure 
of low-energy spin states [28[ . Coherence times of up to 
~ 3 /is 129] which can persist up to relatively high tem- 
peratures of the order of few Kelvin are sensitive to the 
isotopic composition of the molecule. A universal set of 
quantum gates can be applied in a system of coupled an- 
tiferromagnetic ring molecules, without the need for local 
manipulation [301 ] . The presence of many magnetic cen- 
ters with the coupled spins allows for the construction 
of spin cluster qubits that can be manipulated by rela- 
tively simple means [3lj |. In polyoxometalates, the spin 
structure of the molecule is sensitive to the addition of 
charge, and controlled delivery and removal of charges 
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via an STM tip can produce useful quantum gates [32j . 
Chemical bonds between the molecules can be engineered 
to produce the permanent coupling between the molec- 
ular spins and allow for interaction between the qubits 

Sensitivity of molecular state to the addition of 
charge was demonstrated in the tunneling through sin- 
gle molecules [Hj], and used to control the spin state of 
a MN [3^ . Transport studies of the MNs can provide a 
sensitive probe of their spin structure [27l . |36| - [33 | . 

The most straightforward and traditional way of con- 
trolling magnetic molecules is by applying an external 
magnetic field. With carefully crafted ESR pulses, it is 
possible to perform the Grover algorithm, or use the low- 
energy sector of the molecular nanomagnet as a dense 
classical memory [|J. Unfortunately the approaches 
based on magnetic fields face a significant drawback in 
the large-scale quantum control application. Typically, 
the quantum manipulation has to be performed on the 
very short spatial and temporal scales, while the local 
application of rapidly varying magnetic field presents a 
challenging experimental problem. For that reason, the 
schemes for quantum computing tend to rely on modify- 
ing the spin dynamics that is caused by intramolecular 
interaction, rather than on the direct manipulation of 
spins 40]. 

For the applications that require quantum control, the 
electric fields offer an attractive alternative for spin ma- 
nipulation in the molecular nanomagnets [4lj . One major 
advantage is that they can be applied to a very small 
volume via an STM tip [H, [43|, and rapidly turned 
on and off by applying voltage pulses to the electrodes 
placed close to the molecules that are being manipu- 
lated. Switchable coupling between different nanomag- 
nets is essential for qubit implementation. At present, 
this can be implemented only locally, and the interac- 
tion is practically untunable. The use of microwave cav- 
ities can offer a solution to this problem. By placing 
the nanomagnets inside a microwave cavity, one can ob- 
tain a fully controllable, long-range interaction between 
them[4l|. This coupling relies on the presence of a quan- 
tum electric field inside such a cavity, which mediates 
the interaction between distant nanomagnets. The in- 
teraction can be tuned by tuning each molecule in- or 
out-of-resonance with the cavity field using local electric 
or magnetic fields [4lj. The spins, however, do not cou- 
ple directly to the electric fields, classical or quantum, 
and therefore any electric spin manipulation is indirect, 
and involves the modification of molecular orbitals or the 
spin-orbit interaction. 

The description of the molecular nanomagnets in terms 
of spins is an effective low-energy theory that does not 
carry information about the orbital states. However, it 
is still possible to predict the form of spin-electric cou- 
pling from symmetry considerations and single out the 
molecules in which such a coupling is possible. In par- 
ticular, the molecules with the triangular arrangement of 
antiferromagnetically coupled spin- 1/2 magnetic centers 



interact with external electric field through chirality of 
their spin structure [4l|, |44| • The same coupling of chi- 
rality to the external electric field was derived for the 
triangular Mott insulators (45j . 

While the symmetry of a molecule sets the form of 
spin-electric coupling, no symmetry analysis can predict 
the size of the corresponding coupling constant. The 
coupling strength will depend on the underlying mech- 
anism that correlates the spin and orbital states, and on 
the detailed structure of low-energy molecular orbitals. 
To identify molecules that can be efficiently manipulated 
by electric fields, it is necessary to perform an extensive 
search among the molecules with the right symmetries 
and look for the ones that also have a large coupling 
constant. Unfortunately, this search has to proceed by 
ab-initio calculations of the coupling constants for a class 
of molecules of a given symmetry, or by an indiscriminate 
experimental scanning of all of the available molecules. 

In this paper, we contribute to the search for molecules 
that exhibit strong spin-electric coupling. Based on the 
symmetry analysis, we identify the parameters of the spin 
Hamiltonian that can change in the magnetic field, and 
cause spin-electric coupling. We study the mechanisms 
that lead to this coupling and describe the experiments 
that can detect it. 

We will consider the spin electric coupling in the lan- 
guage of effective model, namely either the spin Hamil- 
tonian, or the Hubbard model. In reality the mechanism 
behind the spin-electric coupling involves either the mod- 
ification of the electronic orbitals in an external field and 
the Coulomb repulsion of electrons, or the much weaker 
direct spin-orbit coupling to the external fields. A deriva- 
tion of spin-electric coupling from this realistic picture 
would require the knowledge of electronic orbitals from 
an ab-initio calculation, and the distribution of electric 
field within the molecule. Both of these problems require 
substantial computational power, and can not be per- 
formed routinely. Since the electric field acts primarily 
on the orbital degrees of freedom, and the spin Hamilto- 
nian carries no information about the orbital states, we 
provide a description in terms of a Hubbard model that 
still contains some information about the orbital states. 
We can then described the properties of the molecule that 
allow for strong spin-electric coupling in the language of 
orbitals that offers some intuitive understanding of the 
underlying mechanisms of interaction. 

We identify the response of an MN with spin-electric 
coupling in the standard measurements of ESR, nuclear 
magnetic resonance (NMR), magnetization, polarization, 
linear magnetoelectric effect, and specific heat measure- 
ments. 

In Sec. II we present a symmetry analysis of the spin- 
electric coupling in the ring-shaped molecules with an- 
tifcrromagnetic coupling of spins. In Sec. IV, we de- 
scribe the MNs using the Hubbard model, and relate the 
symmetry-based conclusions to the structure of molecu- 
lar orbitals. In Sec V, we analyze the experimental sig- 
natures of spin-electric coupling, and present our conclu- 
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FIG. 1: (Color online) Schematics of the Si = 1/2 triangu- 
lar molecule in electric field. The antiferromagnetic exchange 
couplings, represented by the bonds with thickness propor- 
tional to Ja+i, are modified in electric field. In the absence 
of electric field, exchange couplings are equal Ju+i = Jjj+i, 
fade colors (grey online). The full color (blue online) triangle 
represents the exchange interaction strengths in electric field. 



sions in Sec. VI. 



II. SYMMETRY ANALYSIS OF 
ANTIFERROMAGNETIC SPIN RINGS 

Spin chains whose ground state multiplet consists of 
two quasi-degenerate S = 1/2 doublets represent suit- 
able candidates for the manipulation of the spin state 
by pulsed electric fields. Such a ground-state multiplet 
characterizes a number of frustrated spin rings, consisting 
of an odd number of half-integer spins. In the following 
we consider prototypical examples of such systems. 



A. Triangle of s = 1/2 spins 

The low-energy properties of most molecular nanomag- 
nets (MNs) are well described in terms of spin degrees of 
freedom alone. Within the spin-Hamiltonian approach, 
the coupling of external electric fields to the molecule can 
be accounted by suitably renormalizing the physical pa- 
rameters. In the following, we use the symmetry of the 
molecules to calculate the changes of spin-Hamiltonian 
parameters, to identify the system's eigenstates, and to 



deduce the allowed transitions. Quantitative estimates 
of the parameters entering the spin Hamiltonian require 
the use of ab-initio calculations 46], or the comparison 
with experiments. The simplest example of a spin system 
which may couple to an external electric field in a non- 
trivial way is a triangle of s = 1/2 spins, like, for example, 
the CU3 MNJStJ. The schematics of such a spin system 
in the presence of an electric field is showed in Fig. [1] 
Its spin Hamiltonian, for the moment in the absence of 
any external fields (magnetic or electric), reads: 

N N 
#spin = ^2 Jii+lSi ■ Si+1 + ^ Djj+1 • (Si X S i+ i), (1) 
i=l i=l 

with N = 3 and S4 = Si in the summation over i. The 
first term in Eq. (fTJ) represents the isotropic Heisen- 
berg exchange Hamiltonian with the exchange couplings 
Ju+i between the spins Sj and Sj+i, and the second term 
represents the Dzyalozhinsky-Moriya (DM) interaction 
due to the presence of spin-orbit interaction (SOI) in the 
molecule, with the DM vectors D,-$+i. The states of the 
spin S — 1/2 triangle can be found by forming the di- 
rect product of the SU (2) representations of three spins 
5=1/2: Dtot = £>(V2)®3 = 2£>( 1 /2) © D^' 2 \ meaning 
there are eight states in total. The point group symme- 
try of the molecule is D3^[13|, i.e. the triangle is assumed 
to be equilateral. The symmetry imposes the fol- 
lowing restrictions on the spin Hamiltonian parameters: 
J li+ i = J and Dfj+j = 0, and Df i+1 = D z . However, 
if lower symmetry is considered these restrictions will be 
relaxed. The spin states in a form adapted to the rota- 
tional symmetry C3 of the system are 

i^ti /2 > = 4?x>Mi+tt) (2) 

Vd 3=0 

1^=3/2) = Ittt), (3) 

where €j = exp (2z7r/3j) and j — 0, 1, 2. The states with 
opposite spin projection M' = —M, i.e. with all spins 
flipped can be written in an identical way (not shown). 
These states are already the symmetry adapted basis 
functions of the point group D 3 /j. Moreover, these are 
eigenstates of the chirality operator 

Cz = ^=si • (s 2 x s 3 ), (4) 

with c z \^J ±l/2 ) - iiv^L/s), ^i<L ±1/2 > = 

and C7 z |'0 M=± 3/ 2 ) = 0. The above states in Eq. (|3|) carry 
different total spin. There are two spin 5 = 1/2 states, 
corresponding to k — 1,2, and a spin 5 = 3/2 state 
corresponding to k — 0. Obviously, the states ^^=±3/2) 
have S = 3/2. 

In an even-spin system, double valued point groups, in- 
stead of single valued groups, are usually used in order to 
describe the states, the splittings and the allowed transi- 
tions (magnetic or electric) [48[ . In the presence of of spin- 
orbit interaction the splittings can be accounted for either 
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by single group analysis (perturbatively) , or by double 
group analysis (exact). In the following, we analyze the 
spectrum and the allowed transitions by both single val- 
ued point group analysis and double valued point group 
analysis. 



E and the bond s± — s 2 - For Cm = (C x , C y ,0) we have 



c- = EO^XV'mI + I^X^mI). ( 8 ) 

M 

Cy = iJ2 (\^mWS\ - I^X^l) • (9) 



1. Single valued group analysis of the s = 1/2 spin triangle 

In the single valued point group T>zh, the states 
I^mL±i/2) w ^ tn ^ = 1j2 form the basis of the two di- 
mensional irreducible representation E , while the states 
; 2 }j and the \4>m=±3/2) transform as A 2 . The al- 
lowed electric transitions in the system are determined 
by the transformation properties of the basis states. 

The simplest and possibly the dominant dependence of 
the spin Hamiltonian on the applied electric field comes 
via the modification of the exchange interactions, like 
depicted in Fig. [TJ This gives rise to the following term 
in the spin Hamiltonian 



SH (E) 



5Ju+i(E) Si ■ Sj+i, 



(5) 



where 5Ju+i (E) w d^ + i -E, with djj+i being vectors that 
describe the electric-dipole coupling of the bond Sj — s i+ i 
to the electric field E in leading order. There are three 
such vector parameters and thus nine scalar parameters 
in total. However, symmetry will allow to drastically 
reduce the number of free parameters by providing rela- 
tions between them. The S — 3/2 states of the unper- 
turbed spin Hamiltonian form the multiplet 4 ^4 2 , while 
the S =1/2 states form two multiplets 2 E . The electric 
dipole Hamiltonian is i? e _d = ~ e J2i^ ' F i = — eE • R, 
with e standing for the electron charge, being the coor- 
dinates of the i-th electron and R = rj. The non-zero 
electric dipole matrix elements of -ff c _d in the D^h sym- 
metric molecule are 



(1,2) 
M 



-ex\tp^'P) 



i(ip 



(1,2) I 
M 



■ey\ip 



(2,1)\ 
M' I 



d5 



MM' 



proportional to the effective electric dipole parameter d. 
The value of d is not determined by symmetry, and has 
to be found by some other means (ab-initio, Hubbard 
modeling, experiments, etc). We mention that all the 
other matrix elements are zero, e.g. 

(12) (12) 

iypM I — e VWM> ) = 0' e ^ c - We see that the electric 
field acts only in the low-energy sector, which allows us 
to write the effective spin-electric coupling Hamiltonian 
acting in the lowest quadruplet as 



o-cff 
H c-d 



dE' ■ Cii 



(7) 



where E' = K z (7n/6 - 26)E, with K z (<f>) describing the 
rotation with an angle <fi about the z axis, and 9 is the 
angle between in-plane component En of the electric field 



The low-energy spectrum in the presence of electric field 
and the related states can be expressed in terms of the 
spin Hamiltonian Eq. ([5]), so that we find anisotropic 
variations of the exchange coupling constants: 



5J ii+ i(E) = — |E|||cosf —1 



(10) 



which depend on the angle 9 and the projection of the 
electric field E on the plane of the triangle. In the Si — 
1/2 triangle the C-operators can be written as 

2 

C x = --(si • s 2 - 2s 2 • S3 + 83 -si), (11) 
2 

C y = -p(Sl • s 2 - s 3 • si), (12) 
V3 

with [Cj , Cj = 2ietjkCk {^ijk are the Levi-Civita 
symbols) (4l[ l45j . From the above relations we can con- 
clude that (i) only the electric field component per- 
pendicular to the bond and lying in the plane of the 
molecule gives rise to spin-electric coupling and (ii) 
there is only one free parameter d describing the cou- 
pling of the spin system to electric fields and du+i = 
Ad/3 (sin (2mt/3), cos (2in/3), 0), where i = 1,2,3 labels 
the triangle sites and 4 = 1. 

The SOI in a D^h symmetric MN is constrained by 
the transformation properties of the localized orbitals. It 
reads 

Hso — AgoT^S^ + Aso (T E » S- +T E »S+), (13) 

with Tr being tensor operators transforming according to 
the irreducible representation r[48j. The non-zero ma- 
trix elements of this SOI Hamiltonian in the low-energy 



W quadruplet read (V>m^ I-^soIV'm' 2 ^) = ±M\I q 8mm' 



so 



that the SOI takes the following effective form 



Hi 



so 



(14) 



with Aso = ^so an< ^ S z — Y^i s l- An effective SOI 
Hamiltonian is obtained also from the DM SOI Hamil- 
tonian in Eq. (JXJ) - The constraints -D^'i = and 
Df i+1 = D z on the DM vectors due to D 3?i symmetry 
of the molecule, give rise to the same effective SOI in Eq. 



(IT41 . with D z = Ag Q . Thus, as expected, the molecu- 
lar SOI and the DM SOI give rise to the same effective 
SOI Hamiltonian acting in the low energy quadruplet. 
Like in the case of the electric dipole parameter d, find- 
ing Z) z (Ag ) requires more than symmetry, like ab-initio 
methods or experiments. The transverse SOI, with inter- 
action strength Ag~ does not act within the low-energy 
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S*— 1/2: a=l 



S,=-l/2; C z =-l 



B-field c, = j^s t ■ (s 2 x s 3 ) B-field 
transitions I 1 transitions 



~~mm transitions 



S 2 =l/2; C 2 =l 



S 2 =l/2; C 2 =-l 



FIG. 2: The spin transitions in the Sj = 1/2 triangle induced 
by electric and magnetic fields. The electric field causes tran- 
sitions between the states of opposite chiralities Cz and equal 
spin projections S z (horizontal arrows), while the magnetic 
field instead causes transitions between the states of opposite 
spin projections S z and equal chiralities C z (vertical arrows). 



space, and its effect will appear only in higher orders of 
perturbation theory in 1/J. 

An external magnetic field couples to the spin via the 
Zeeman term Hz = B • gS, with g = diagjgy, gy, g±} be- 
ing the g- factor tensor in Dsh. The full effective Hamilto- 
nian describing the low-energy quadruplet in the presence 
of SOI, electric field and magnetic field read 



H cS = A so C z S z + B • gS + dE' • C,, 



(15) 



Note that [C, S] = 0, and chirality and spin act as in- 
dependent spin 1/2 degrees of freedom. Furthermore, in 
the absence of SOI the chirality C and the spin S evolve 
independently. However, the SOI couples the two and 
provides with means for electric control of both spin and 
chirality. Vice-versa, magnetic fields can also couple to 
chirality due to SOI. Also, while magnetic fields (time- 
dependent) cause transitions between states of opposite 
spin projection M but with the same chirality C z , the 
electric field does the opposite: it causes transitions be- 
tween states of opposite chirality C Zl but carrying the 
same M. Full control of the lowest quadruplet is thus re- 
alized in the presence of both electric and magnetic fields, 
as can be seen in Fig. [2] 



2. Double valued group states of the s = 1/2 spin triangle 

The double group representations allow to non- 
perturbatively describe the magnetic and electric transi- 
tions in the presence of spin-orbit interaction. The lowest 
quadruplet consists of two Kramers doublets, one of them 
transforming like E ~ ( | — 1/2) , 1 1/2) ) , and the other one 
according to E" ~ (| - 3/2), |3/2)). Here (|M), | - M)) 



represent pairs of eigenstates of a given angular momen- 
tum J > M, with spin projection ±M. For example, 
if M = 1/2, then J = 1/2, 3/2, .... The higher energy 
states instead (S = 3/2 states), transform now not as 
A' 2 , but as E (M = ±1/2) and as E" (M = ±3/2). 
Thus, the 5 = 1/2 states mix with the 5* = 3/2 states, 
but only the ones transforming according to the same 
representations, i.e. there is no mixing between E and 
E due to spin-orbit interaction. The magnetic dipole 
transitions take place between E and E , and within 
E and E , respectively, while electric dipole transitions 
take place only between E and E . The selection rules 
for the electric transitions are AM = ±2, while for the 
magnetic transitions these are AAf = ± 1 . We see that 
within the double group analysis, i.e. in the presence 
of SOI, there are allowed electric dipole transitions also 
within the 5 = 3/2 subspace. 

Using both the single group and double group analy- 
sis we can pinpoint to the transitions that arise in the 
absence or only in the presence of SOI. Therefore, the 
electric dipole transitions present in the single-group are 
a consequence of the modified exchange interaction, and 
can arise even in the absence of SOI, while the ones that 
show up only in the double group analysis are a con- 
sequence of the SOI (or modification of SOI in electric 
field). 

We now can establish several selection rules for the 
SOI, electric field and magnetic field induced transitions. 
Note that the above analysis was exact in SOI. However, 
it instructive to treat electric field, magnetic fields and 
SOI on the same footing. First, we find that the electric 
dipole transitions fulfill the selection rules AC Z = ±1 and 
AS Z = 0, meaning that electric field only couples states 
within the lowest quadruplet. The SOI transitions show 
a richer structure. We can separate the SOI interaction 
in two parts: the perpendicular SOI, quantified by D z in 
the DM interaction Hamiltonian, and the in-plane SOI, 
quantified by D x ^ y in the DM interaction Hamiltonian, 
respectively. By doing so, we find that the D z SOI terms 
obey the selections rules AC Z = and AS Z = 0, while 
for the D XtV terms we get the selection rules AC Z = ±1 
and AS Z = ±1. We see in-plane SOI {D x<y terms) do 
not cause any splitting in the ground state and can lead 
to observable effects only in second order in perturba- 
tion theory in D x y / J . Also, note that if an symmetry is 
present, D x ^ v = and thus there are no in-plane SOI ef- 
fects at all. Modification of these terms due to an in-plane 
external electric field E, however, lead to different selec- 
tion rules: changes of D z terms lead to AC Z — ±1 and 
AS Z = 0, while modification of D x ^ y lead to AC 2 = 0, ±2 
and AS Z = ±1. The magnetic field transitions obey 
the selection rules AS Z — 0,±1 and AC Z = 0. Thus, 
we can make clear distinction between pure electric field 
transitions, SOI-mediated electric transitions and mag- 
netic transitions. This distinction between the electric 
and magnetic field induced transitions could be used to 
extract the spin-electric coupling strength parameter d 
from spectroscopic measurements. 
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B. Spin s = 3/2 triangle 

The spin s = 3/2 triangle has a more complex level 
structure than the s = 1/2 triangle due to its higher 
spin. The spin Hamiltonian, however, is similar to the 
one in Eq. JTJ for s = 1/2, and the reduction of the rep- 
resentation of three spins S = 3/2 is D to t = £>( 3 / 2 )® 3 = 
2 £>(i/2) © 4i^( 3 / 2 ) © 3D^ © 2£>( 7 / 2 ) © £( 9 / 2 \ a total 
of 64 spin states. The total number of irreducible rep- 
resentations is the same as in the s = 1/2 case, and we 
need only to identify these basis states in terms of the 
spin states. The s = 3/2 triangle states can be defined 
according to their transformation properties under three- 
fold rotations C3 in D^h and are of the following form 

= P fe 3 |M,i), (16) 

Pi - TTfE^, (17) 

3=0 

where = exp (2i7rj'fe/3), C| are the 3-fold rotation of 
order j, and j, k = 0, 1, 2. The states \M, i) = \a1a2cr3) 
represent all possible states (i states in total) with a given 
spin projection M (= <jfc) that cannot be transformed 
into each other by application of the rotation operator 
C3 . These states are showed in Table fl] 



\ i 
M\ 


1 


2 


3 


4 


1/2 


lltt) 


lftw> 


\m) 




3/2 


IMft) 


\xn) 


\m) 


Ittt) 


5/2 


itrt) 


\m) 








7/2 


nu) 











9/2 


Iftftft) 












TABLE I: Non-symmetry adapted states of the s = 3/2 spin 
triangle. We use | ft (JJ.)) = | ± 3/2). 

The corresponding states with all spins flipped, namely 
with M = —M, can be written in a similar form (not 
shown). Having identified the symmetric states in terms 
of the spin states, we proceed to analyze the allowed tran- 
sitions induced in the spin systems by magnetic and elec- 
tric field, both within the single valued group and double 
valued group representations. 

1. Single valued group states of the s = 3/2 triangle 

The above states are basis of the point group D 3 /j, 
but not eigenstates of the total spin operator S 2 , i.e. 
they do not have definite total spin. However, linear 
combinations of states of a given total spin projection 
M and a given 'chiral' numbers k become eigenstates 
of S 2 . The total spin eigenstates can be written as 

IV^m) = J2i(M) a fc,il^M )> where K M ) is the number of 
different states with a given M. The coefficients ak,i are 



to be identified so that these states satisfy S 2 \ip s M ) = 

5(5 + l)\ip { s k) M ), with 5 = 1/2,3/2,5/2,7/2,9/2. The 
states with k — are all transforming according to the 
A 2 representation, while the states with k = 1, 2 are or- 
ganized in doublets, being the bases of the two dimen- 
sional representation E . However, as mentioned above, 
different combinations of symmetry adapted states carry 
different total spin 5. The magnetic and electric transi- 
tions are similar to the ones in the s = 1/2 triangle, in 
the absence of SOI. The electric field causes transitions 
only between states with the same M and 5, but oppo- 
site chirality C z — 7^7j s i 1 ( s 2 x s 3) (this is different from 
the triangle with Sj = 1/2 spins in each of the vertices). 
As for the s = 1/2 spin triangle, there are electric dipole 
transitions within the spin system even in the absence of 
SOI. The ground states is four-fold degenerate consisting 
of two 5=1/2 eigenstates 

- (ei-e 2 )(|^! ) 1/2 >-|^'l ) 1/2 >)),(18) 

+ (^i-e 2 )(|^ V2 )-|^f 1/2 >)).(19) 

We see that, as opposed to the s = 1/2 triangle, the 
lowest states are given by linear combinations of the sev- 
eral M = 1/2 symmetry adapted states (the M = —1/2 
states are obtained by flipping the spins in the states in 
Eqs. (fT5|). (fT9"]). This, however, does not modify the con- 
clusions regarding the electric and magnetic transitions 
in the absence of SOI, these being given by the same rules 
as in the 5=1/2 triangle: electric- field induced transi- 
tions between the states of opposite chirality C z and the 
same spin projection M. The lowest states are still or- 
ganized as spin and chirality eigenstates that are split in 
the presence of SOI as in the previous case. 

In the original spin Hamiltonian in Eq. (flj the electric 
field causes modification of the spin Hamiltonian param- 
eters. As for the spin s = 1/2 triangle, the strongest 
effect comes from modification of the isotropic exchange 
interaction, so that 

3 

«fo(E) = 5>J ii+ i(E) S i-s i+1 , (20) 

i=l 

with 5Ju + i(E) — dE cos (2iri/3 + 9), where is the angle 
between the projection of the external electric field E to 
the molecule's plane and the s\ — s 2 bond, and i = 0, 1, 2. 
The effect of the electric field on the lowest quadruplet 
is found to be similar to the spin s = 1/2 case. While 
the SOI splits the two chiral states without mixing them 
(at least in lowest order), the electric field, on the other 
hand, mixes the chiral states. The effective Hamiltonian 



acting in the lowest quadruplet reads 

H eS = A so C z S z + B-gS + d'E- C,,. (21) 
Above, d' = 3d/2, C\\ = {C x ,C y ,0), with 
C x = Em\^mWm ] \ + \^mWm ] \ and C x = 
i EaM { m } )(^m I " \4>m)(^m I), and A so stands for the 
SO splitting. However, in this situation the in-plane chi- 
rality operators C x , y cannot be written in a simple form 
as a function of the individual spin operators, as opposed 
to the s = 1/2 triangle. 

2. Double valued group states of the s — 3/2 triangle 

The double group representation allows to identify the 
couplings between different spin states induced by the 
SOI and to identify the allowed magnetic dipole transi- 
tions. Due to SOI, the electric field induced spin transi- 
tions will take place also outside the spin quadruplet. In 
the absence of extra degeneracies (induced, for example, 
by external magnetic fields), however, these transitions 
are strongly reduced due the gap of the order J. We 
can then focus, as for the 5=1/2 triangle, only on the 
lowest quadruplet. These states are organized in two 
Kramer doublets of the form (|M), | — M)), one trans- 
forming as E ~ (| 1/2), | — 1/2)) and the other one as 
E" ~ (| - 3/2), |3/2)). Here again, (\M), \ - M)) rep- 
resent angular momentum J < M eigenstates with spin 
projection ±M. 

As in the case of the s = 1/2 triangle, the electric 
field induced transitions take place between E and E , 
with the selection rules AM = ±2. Magnetic transitions 
instead take place both within and between E and E , 
satisfying the selection rules AM = 0, ±1. 

If we now treat the SOI, electric field and magnetic 
fields on the same footing, we arrive at the same selection 
rules as for the s = 1/2 triangle, namely AC Z — ±1 and 
AS Z = for electric transitions, AC Z = 0, ±1 and AS Z = 
0, ±1 for SOI transitions, and AC Z = and AS Z = 0, ±1 
for magnetic transitions, respectively. 

C. Spin s = 1/2 pentagon 

We now analyze the spin-electric coupling in a pentag- 
onal molecule with a spin s = 1/2 in each of the vertices, 
like depicted schematically in Fig. |3J As in the case of 
the spin triangle, an external electric field E gives rise 
to modification of exchange interaction Ju+i in Eq. fT]). 
However, the net spin-electric coupling in the lowest spin 
sector can only be mediated by SOI. i.e. via the DM in- 
teraction (which can be also modified in the presence of 
the E-field). 

To make the analysis simpler, we assume in the fol- 
lowing that the pentagonal spin molecule possesses a D 5 
point group symmetry, thus no horizontal reflection plane 
ah- However, no generality is lost, since lower symmetry 
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FIG. 3: (Color online) Schematics of a pentagonal spin ring 
molecule in electric field E, light (green) arrow. The molecule 
in the absence of electric field is depicted in fade colors, while 
the full colors represent the molecules in electric field. Thick- 
ness of the bonds represents the strength of antiferromagnetic 
exchange interaction between the spins. An electric field mod- 
ifies the strengths of spin exchange couplings J«+i. 



implies more allowed transitions in the spin system. If, 
for example, in the lower symmetric situation some tran- 
sitions are forbidden, these transitions will be forbidden 
in the higher symmetry case. The Hamiltonian is given 
in Eq. ([TJ with N = 5. The states of the pentagon 
are found from the product of the individual spin repre- 
sentations D tot = D^® 5 = 5D^ © AD^ © D^' 2 \ 
meaning there are 32 spin states in total. As before, these 
states can be organized in a symmetry adapted basis in 
the following way 

= P*\M,i), (22) 
P " = 7fE4' C i> (23) 

where e| = exp (2iirjk/5) with k,j = 0, ...,5, are 
the 5-fold rotations of order j. The states |M, i) = 
| cr i <72 03 a 4 05) represent all possible states (i states in 
total) with a given spin projection M(= X)fe cr fe) that 
cannot be transformed into each other by application 
of the rotation operator C|. These states are showed 
in in Table |TT] and the corresponding states with all 
spins flipped, i.e. M —> —M states (not shown). In 
the absence of SOI there is no mixing of different k 



8 



\ i 
M\ 


1 


2 


1/2 


mm) 


imtt) 


3/2 


utttt) 





5/2 


ittttt) 






TABLE II: Spin s = 1/2 pentagon non-symmetry adapted 
states. 



that d is in principle non zero in D5 point group sym- 
metry. However, the matrix elements of the spin-electric 
Hamiltonian within the lowest quadruplet are all zero, 
i.e. {i} { gl 1/2 M \5H c _ d {^)\^J 1/2 M ,) = 0. This means 
that electric field has no effect on the lowest quadru- 
plet, as found out also by purely symmetry arguments. 
Therefore, we may expect that the spin-electric coupling 
in pentagonal spin molecule is caused by SO effects. 



states, i.e. the chirality is a good quantum number. 
In this case the chirality is quantified by the operator 

C z = 1/(2V5 + 2y/E)J2i s i ■ ( s i+i x s »+2) (the prefac- 
tor is chosen for convenience; see below). As in the 
s = 1/2,3/2 spin triangles, the above states are not 
yet the eigenstates of the Hamiltonian and we have 
to solve the equation S 2 \ip^) = S(S + 1)1^), with 

1^4 ) = Efc(M) a k,i\^M >• The ground state is spanned, 
again, by four states, two Kramers doublets with spin 
5 = 1/2. In the following we inspect the level struc- 
ture of these four states in terms of the above symmetry 
adapted states. 

1. Single valued group s — 1/2 pentagon 

We focus here only on the four lowest energy states, 
which are two pairs of S — 1/2 states. The first (second) 
pair is given by linear combination of states with chirality 
k = 1 (k = 4) and spin projection M = ±1/2. We obtain 

\vb {k) ) - 1 ( 1 k- (M) ) 

I^S=l/2,M=±l/2/ — /5I /2fc7r\ M=±1 / 2 

2cos {— ) 

+ 2e §c OS (^)|^ 1/2 >), ( 24 ) 

so that C z \tp^ =±1 , 2 ) = (— 1)' c |V'mL±i/9)- These states 
(for a given M projection) form the basis of the two di- 
mensional irreducible representation E±. We are now in 
positions to investigate the allowed electric dipole tran- 
sitions within this lowest subspace. The in-plane elec- 
tric dipole d = (d x ,d y ) forms a basis of the irreducible 
representation E\ in D5. By calculating the product 
Ei®Ei®E\ = 2E\®2E 2 we see that the totally symmet- 
ric representation A\ of D5 is absent. Therefore, there 
are no electric dipole transitions within the four dimen- 
sional subspace in the absence of SOI. 

As in the previous two cases, the coupling of the spin 
Hamiltonian to electric field comes via modification of 
the spin Hamiltonian parameters. If only the modifica- 
tion of the isotropic exchange Hamiltonian is taken into 
account, the spin-electric Hamiltonian takes the same 
form as in Eq. 0, with 5J ii+1 (E) = dE cos (2«7r/5 + 0), 
i = 1 . . . 5. The parameter d quantifies the electric dipole 
coupling of each of the bonds and 9 is the angle be- 
tween the electric field E and the bond Si — S2. Note 



2. Double valued group s = 1/2 pentagon 



Double valued group analysis allows identifying of the 
level structure and the allowed transitions in the pres- 
ence of SOI and magnetic fields. The lowest four states 
in the double group D 5 are described by the two dimen- 
sional irreducible representations E 1 ~ (| — 1/2), |l/2)) 
and E 1 ~ (| — 3/2), 1 3/2) ) , respectively. Since both 
the magnetic /z and electric d dipoles transform as E\ 
in D 5 , both electric and electric transitions will take 
place between the same pair of states. The products of 
the irreducible representations that labels the states in 
the low-energy quadruplet read: E 1 ® E 2 = E\ ® E2, 
E[ <g> E[ = A x © A 2 © E 1 and E 2 <g> E 2 = A x © A 2 © E 2 . 
These equalities imply the same selection rules in the 
lowest subspace as for the spin triangle case: AM = ±2 
(| ± 1/2) O I =F 3/2) ) for electric dipole transitions, and 
AM = ±1 (| ±1/2) <-> I ±1/2) and | ± 1/2) O | =b 3/2) ), 
for the magnetic ones. 

The main feature of pentagonal spin ring is the ab- 
sence of electric dipole transitions in the lowest quadru- 
plet in the absence of SOI. This is to be contrasted to 
the spin triangle case, where spin-electric coupling ex- 
ists in the ground state even in the absence of SOI. This 
feature finds its explanation from the interplay between 
the selection rules for electric field transitions and the 
ones for the SOI. In fact, these selection rules are by 
no means different from the triangular spin rings. Since 
the ground state is spanned by four states with chirality 
C z = 1,4 and spin S z — ±1/2, we see that the condi- 
tion AC Z = ±1 for the electric field transitions implies 
no electric field coupling within the ground state! In 
the presence of SOI though, spin electric coupling is still 
possible, but it will be (D x ^ y / J) times smaller than in 
triangles. Spin-electric coupling can arise also via modifi- 
cation of the DM vectors D x ^ ViZ in electric field. However, 
the selection rules for this transitions are, like for the tri- 
angle, AC Z — 0, ±2 and AS Z = 0, ±1. This means direct 
splitting in the ground state, and thus we expect that 
for pentagon spin ring the electric dipole response will 
be much weaker. 
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III. HUBBARD MODEL OF A MOLECULAR 
NANOMAGNET 

Spin-Hamiltonian models of molecular nanomagnets 
are based on the assumption that the spins on magnetic 
centers are the only relevant degrees of freedom. This as- 
sumption of fully quenched and localized orbitals allows 
for the relatively simple predictions of spin structure in 
the low-energy states of the molecule. However, since 
the orbital dynamics plays a crucial role in spin-electric 
coupling, spin-Hamiltonian models are unable to predict 
the corresponding coupling constants. In this Section, we 
relax the assumption of quenched and localized orbitals 
and treat the orbital degrees of freedom of electrons on 
magnetic ions within a Hubbard model. This provides 
an intuitive picture of spin-electric coupling in terms of 
the deformation of the molecular orbitals induced by the 
external field. Besides, in the limit of strong quench- 
ing of the orbitals, the Hubbard model reproduces a spin 
Hamiltonian, similar to the results found in the studies 
of cuprates [49Tj5l| and multiferroics [H, H3| ■ I* 1 partic- 
ular, we find the relation between modifications of the 
electronic hopping matrix elements induced by the field 
and that of the spin-electric coupling in the spin Hamil- 
tonian, thus providing a guide for the estimate of the size 
of spin-electric coupling in a molecule. 

The outline of the present Section is the following. In 
Subsection IIII A[ we introduce the Hubbard model of a 
spin chain with the shape of regular n-tangon, and de- 
rive the resulting symmetry constraints for the hopping 
parameters. In Subsection IIII Bl we assume a direct elec- 
tron hopping between magnetic sites, and derive the spin 
Hamiltonian of a spin triangle from the Hubbard model, 
in the limit of large on-site repulsions; we thus express 
the coupling to electric fields in terms of the Hubbard- 
model parameters. In Subsection IIII CI we introduce a 
Hubbard model of a magnetic coupling in the case where 
this is mediated by a non-magnetic bridge between the 
magnetic centers; also in this case, we find a connection 
between the modification of the bridge and spin-electric 
coupling. 



A. Parameters of the Hubbard model of molecular 
nanomagnets 

Magnetic properties of molecular nanomagnets are 
governed by the spin state of few electrons in the highest 
partially occupied atomic orbitals, split by the molecu- 
lar field. The spin density is localized on the magnetic 
centers [54[ , and thus the low-energy magnetic properties 
are correctly described by quantum models of interacting 
localized spins [H, [56[ . 

The response of molecular nanomagnets to electric 
fields, as a matter of principle, does not have to be gov- 
erned by the electrons occupying the same orbitals that 
determine the molecule's spin. However, the quantum 
control of single molecule magnets by electric fields de- 



pends on the electrons that both react to electric fields 
and produce the magnetic response. Therefore, the mod- 
els of molecular nanomagnets that consider only few or- 
bitals can provide useful information about the electric 
control of spins. 

Hubbard model provides a simplified description of or- 
bital degrees of freedom by including only one or few 
localized orbitals on each magnetic center. Furthermore, 
the interaction between electrons is accounted for only 
by introducing the energies of the atomic configurations 
with different occupation numbers. The Hubbard model 
of the MN is given by: 



J2 Yl c i« (*<W + ^ • OYtf ) c jp + h. c. 



(25) 

where cj CT (c J(T ) creates (annihilates) an electron with 
spin a =T, J, on the orbital localized on jth atom, and 



riju = c'j a Cj a is the corresponding number operator. 



Model parameters Uj, describe the energy of spin 
up (down) electrons electrons on the site j. Hopping pa- 
rameters Uj , Py describe the spin- independent and spin- 
dependent hopping between sites i and j. 

We assume that the largest energy scale is the split- 
ting between the energy of the highest occupied atomic 
orbital and lowest unoccupied one, induced by the molec- 
ular crystal field: this justifies the inclusion of one orbital 
only for each magnetic center. The on-site repulsion en- 
ergy is the next largest energy scale in the problem, be- 
ing Uj larger than the hopping coefficients. Amongst 
these, processes involving states of different spin, medi- 
ated by spin-orbit interaction, are described by the x and 
y components of Py. The parameters Pij- Z , instead, de- 
scribe the difference of the hopping matrix elements be- 
tween spin-up and spin-down electrons. In the following, 
we shall consider both the case where electron hopping 
takes place directly between neighboring magnetic ions 
and that where the magnetic interaction is mediated by 
bridges of non-magnetic atoms. The Hubbard Hamilto- 
nian can be approximated by a spin Hamiltonian model 
in the limit |£y|, |P»j| *C Uj. The symmetry constraints 
on the spin Hamiltonian parameters can be deduced from 
those on the Hubbard model parameters [49| . If the spin- 
independent hopping dominates (\t\ » |P|), the resulting 
spin Hamiltonian will contain the Heisenberg exchange 
terms and a small additional spin-anisotropic interaction. 
If \t\ > |P|, the size of spin-dependent interactions in the 
spin Hamiltonian will be comparable to the Heisenberg 
terms. Both these cases appear in the molecule nano- 
magnets [3 53 13, i| . 

Symmetry of the molecule imposes constraints to the 
Hubbard model, thus reducing the number of free pa- 
rameters. The on-site repulsion parameters Uj are equal 
for all equivalent magnetic ions. In the molecules of the 
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form of regular n-tagon, all of the spin- independent hop- 
ping parameters are equal, due to the C n symmetry. The 
spin-dependent hopping elements are related by both the 
full symmetry of the molecule and the local symmetry of 
localized orbitals. For example, in the case of localized 
orbitals in a regular polygon that are invariant under the 
local symmetry group of the magnetic center, 



P, 



exp 



.2n(j-k) 



Pk,k+l;x, 



(26) 



with the convention that site n + 1 coincides with site 
1. In this case, there is only one free parameter that 
determines all of the P x matrix elements. Therefore, 
the regular n— tagon molecule in the absence of exter- 
nal electric and magnetic fields can be described by a 
Hubbard model, with five independent parameters: U, t, 
P12. In addition, the a v symmetry, if present will impose 
P12 = pe z , thus reducing the number of free parameters 
to three. 



with <jj — a for j ^ k and aj = a, for j = k, and (ii) the 
double-occupied sites 



(30) 



with k = 1, 2, 3 and p ^= k. 

The states in Eqs. (|28|) . (|35|) and (f30|) are degenerate 
with energies E = eo, E — 3eo and E = 3eo + U , re- 
spectively. Note that these state are eigenstates of the 
Hamiltonian in Eq. (f27|) only in the absence of tunneling 
and SOI. 

The above defined states are not yet adapted to the 
symmetry of the system, i.e. they are not basis states 
of the corresponding irreducible representations of T)j,h 
point group. Finding these states is required by the 
fact that the symmetry of the molecule is made visi- 
ble through the hopping and SOI terms in the Hubbard 
Hamiltonian. This is accomplished by using projector 
operators [4^1 . We obtain for the one-electron symmetry 
adapted states. 



B. Hubbard model of the spin triangle: direct 
exchange 

In this Subsection we give a brief description of the 
Hubbard model for a triangular molecule with sym- 
metry. In this model we assume only direct coupling be- 
tween the magnetic centers, thus no bridge in-between. 
Even so, this simplified model catches the main features 
of the effective spin Hamiltonian and gives the micro- 
scopic mechanisms for the spin-electric coupling. The 
Hamiltonian describing the electrons in the triangular 
molecule reads 



H H = 



E J 



c] a (t + ia\ S o)ci+i,a + h.c. 



(27) 



where Aso = P 



e z is the spin-orbit parameter 



(only one), eo is the on-site orbital energy, and U is the 
on-site Coulomb repulsion energy. As stated before, typ- 
ically Asoj \t\ <§C U, which allows for a perturbative treat- 
ment of the hopping and spin-orbit Hamiltonians. These 
assumptions agree well with the numerical calculations 
performed in [54j . 

The perturbation theory program involves the unper- 
turbed states of the system. The first set of unperturbed 
states are the one-electron states 



(28) 



while the three-electron states split in two categories: (i) 
the site singly occupied states 



i^>=n c k-i°>' 



(29) 



^E^io, 

v i—i 



(31) 

(32) 
(33) 



where A' 2 and E' ± are one-dimensional and two- 
dimensional irreducible representations in T>^) ll respec- 
tively. Similarly, the symmetry adapted states with the 
singly-occupied magnetic centers read: 



71 

V 1=1 



(34) 
(35) 



while the symmetry adapted states of the doubly- 
occupied magnetic centers read: 



2 a 
A' . 



^X>&>±w 3 >), 



(36) 



KV = 4E472 1 (I^) + I^}), (37) 

* i— 1 

= 7i E A^m) - (38) 



The tunneling and SOI mixes the singly-occupied and 
doubly-occupied states. Since both the tunneling and 
SOI terms in the Hubbard Hamiltonian transform as 
the totally symmetric irreducible representation A! x in 
D3/j, only states transforming according to the same irre- 
ducible representations T mix. We obtain the perturbed 
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FIG. 4: Spin-Hamiltonian limit. Expectation values of chi- 
rality (C z ) (full lines) and the their bounds of uncertainty 
(C z ) ± AC Z (dotted lines), see text, in the low-energy states 
of the Hubbard model, as a function of the on-site repulsion 
U, at the fixed hopping matrix element t — 1 (left scale). 
The dashed line shows dependence of the double occupancy 
probability in the ground state on the right scale. The spin- 
Hamiltonian description becomes accurate in the U — > oo 
limit. The approach to this limit is slow, and the double 
occupancy probability is proportional to t/U. 



trie field. The second mechanism is due to modification 
of the hopping parameters tu+i in electric field and gives 



Z t *H 

I. a 



Ci+lc 



(42) 



where tf i+1 = ($^1 — er • E^^i^) are new hopping pa- 
rameters induced solely by the electric field E, and $, CT 
are the Wannier states localized on the magnetic cen- 
ters. We can write the E-induced hoppings as = 

Y^ q =x,y,z 1u+i E q> witn 1a+i = ~ eq\® i+la ) being 

electric dipole matrix elements between the i and i + 1 
ions. These matrix elements are not all independent, 
symmetry alone reducing drastically the number of in- 
dependent electric dipole parameters. In order to find 
suitable independent free parameters, we switch from 
the description in terms of localized Wannier orbitals 
<I>i CT , to the description in terms of symmetry adapted 
states , namely from qu+i to q rT > = (4>ra-\q\4>r' a) > where 
r = A 1 ,E ± . In the basis of symmetry adapted states, 
the components q rT > satisfy a number of relations. In 
particular, we find: 



in first order in t/U and Aso : 



(39) (0 



(e-l)(i±ffA S o) 



V2U 



3e(t±a\ so ) 
V2U 



2a 
E'r 



(40) 



Doubly occupied states become high in energy when 
\t\/U, Xso/U <C 1. In this limit, the orbital states are 
quenched into singly-occupied localized atomic orbitals, 
and low-energy behavior is determined by spin and de- 
scribed by a spin Hamiltonian. In this limit the states in 
Eq. (I55|) are exactly the same chiral states in the spin 



Hamiltonian, i.e \ip 



EL 



(1,2) \ 



and 



(0)\ 



The probability of finding two electrons at the same site 
decays as 1/U. The lowest energy states have total spin 
S = 1/2 and the chirality C z = ±1, and the fluctua- 
tions of chirality AC Z = y/{Cf) - (C z ) 2 in the eigen- 
states vanish, see Fig. 21 The chiral states emerge as the 
eigenstates in the large- U limit, when the system is well 
described by the spin Hamiltonian. 

The coupling of the molecule to an external electric 
field E takes place via two mechanisms. The first one 
implies modification of the on-site single particle ener- 
gies eo and leads to the following electric-dipole coupling 
Hamiltonian 



H 



-d 



EE y a | 



Cle 



E V 

V3 



E, 



J 

2<T 



C2o 



E r 



E y 
V3 



c 3<x C 3o 



(41) 



A ,| - ex\ 

E> I - ex V. 
E'J- eX \ l 

ex\i 



'EL, 



(<t>\\- ey\<P\) = (4> E > + \- ™\<t>E' + ) ^«3) 
(4%, + 1 - ey\F E , + ) = (d>%,_ | - ey\4>%-_) 
-i{4>l, + \-ey\P EL )=d EE (45) 

~ ex\cj> EL ) = -i{<t>\\-ey\<t> a E> + ) 
i{<t>% \-evWk' ) = d AE- (46) 



These relations reduce the number of free coupling con- 
stants to two, namely d EE and cIae- 

It is instructive to write first the relation between the 
second quantized operators c i(J (ci a ) and A ra (cr a ), which 
create (annihilate) electrons in localized and symmetry 
adapted states, respectively: 



(A 



la 

t 

'-2a 




(47) 



With these expressions at hand, we can write the electric 
dipole Hamiltonian together with the spin-orbit Hamil- 
tonian in the following form: 



H 



-d 



lZ ( Ec E ' + a C A[a ~ eEA E ,_ a C A[c 



eEc E ,_ a c E , +(7 ) + H.c, (48) 



H c-d = ^d AE {Ec\ {<j c E , +a - Ec\, ia c E ,_ a ) 



Ed EE c% c E > a + H.c. 



(49) 



with a being the geometrical distance between the mag- _ ./q\__ n (A r ^,^ — A A (^(W 

netic ions and E xy the m-plane components ol the elec- 
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where E = E x + iE y (E = E x — iE y ). The symmetry 
adapted states can also be expressed in terms of the sym- 
metry adapted operators c r . The expressions for these 
states are shown in Appendix A. Using these states, we 
can compute all the matrix elements corresponding to the 
electric dipole and SOI Hamiltonian, respectively. The 
explicit form of these matrix elements can be found in 
Appendix B. 

We now compute the electric dipole matrix elements 
between the perturbed chiral states of the E'± sym- 
metry. The question is to what order in t/U and/or 
eEa{d,EEi dAE)/U we want to do it. We use the relations 
\ea\ 3> (Iee, d,AE, which hold in the case of localized or- 
bitals. This leads us to the following matrix element of 
the electric dipole in the ground state: 



Co,, -» C 



\K a ' l^c-dl* 



la 

E' 



— Ed.EE 



(51) 
(52) 



We now relate the SOI matrix elements to the DM vec- 
tors in the effective spin-Hamiltonian. In D^h symmetry, 
the DM term reads 



iD 



i=l 



(53) 



which gives rise to the following non-zero matrix ele- 
ments, 



|Ifso|*& 



± 



V3D Z . 



sign(o-), (54) 



and allows us to make the following identification 



5Asot 
U ' 



(55) 



We see that this SOI term acts exactly as the 'micro- 
scopic' SOI derived before: it splits the chiral states, but 
it does not mix them. 

The Hubbard model with spin-orbit coupling can re- 
produce the energy-level structure of the spin Hamilto- 
nian. In the limit of strong on-site repulsion \t\/U <C 1, 
the atomic orbitals in the triangle vertices are occupied 
by one electron each. The lowest energy manifold con- 
sists of four states with the total spin S^ot = 1/2- These 
states are split from the next four-level Stot = 3/2 man- 
ifold by a gap of the order of t 2 /U . 



C. Superexchange in molecular bonds 

In this Subsection, we use the Hubbard model to de- 
duce the dependence of the spin Hamiltonian of MNs on 
the external electric fields in the case where the coupling 
between magnetic sites is mediated by a non-magnetic 
bridge. In particular, we study how the parameters of 



C 2v —■ C lh 




FIG. 5: Geometry of the bond and reduction of symmetry, 
(a) Electric field E in y direction, leaves the C2v symmetry 
unbroken, (b) An electric field E in z-direction, normal to the 
bond plane, reduces the symmetry to {E, a v }. (c) An electric 
field E in x-direction, along the line connecting the magnetic 
centers, reduces the symmetry to {E,<7h}- (d) In an inho- 
mogeneous staggered electric field E, the reduced symmetry 
group is {E, R v ,tt}- 



the effective spin Hamiltonian depend on the hopping 
matrix elements that are modified by the presence of an 
electric field. This method was successfully applied in 
the studies of strongly correlated electrons, like cuprates 
[Ell ] and multiferroics [53j |. 

In order to describe the magnetic coupling, we consider 
a pair of sites corresponding to the magnetic centers and 
a bridge site. Since the direct overlap of the orbitals lo- 
calized on the magnetic centers is small, we set the direct 
hopping between the magnetic centers to zero, but allow 
for the hopping of electrons between the magnetic sites 
and the bridge site. This hopping gives rise to superex- 
change interaction between the spins on the magnetic 
sites [43|. In the limit of strong on-site repulsions, the 
effective Hamiltonian in the lowest energy sector of the 
bond corresponds to a spin Hamiltonian where the cou- 
pling strengths are determined by the Hubbard model 
parameters. This correspondence provides an intuitive 
picture of the mechanism that leads to the interaction 
between the spins. It also allows us to infer the prop- 
erties of the molecule that lead to a strong spin-electric 
coupling, e.g., the delocalization of the orbitals and their 
local symmetry. 

The Hubbard Hamiltonian of the bond is given by 



i,a0 



tiSa/3 



iPi 



■ Cup I bp + 



h. c. 



(56) 



+ U 1 (n 1 ) + U 2 (n 2 ) + U h (n h ), 



where the indices 1 and 2 refer to the magnetic sites, and 
b refers to the bridge site. We derive the spin Hamil- 
tonian by fourth-order Schrieffer- Wolff transformation of 
the Hamiltonian (jBl)]) . 
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The Schrieffer-Wolf transformation [59j of the bond 
Hamiltonian Hb = Hq + H tun (|56p . where the unper- 
turbed Hamiltonian Hq — Ux(nx) + 1/2(712) + Ub(nb) 
produces an effective low energy Hamiltonian H12 that 
approximately describes the low-energy dynamics of the 
bond. The effective Hamiltonian is 



symmetry are: 



H u = Ve b H h e' b V, 



where the antiunitary operator S is chosen so that the 
low-energy space of Hq is decoupled from the high- 
energy space. This operator is found iteratively, S = 
S^ 1 -* + + . . ., so that the nth order transformation 
g(v) removes the terms that couple the low- and high- 
energy states up to order n. The projector V projects 
to the low-energy states. In our system, the lowest or- 
der Schrieffer- Wolff transformation that gives a nontriv- 
ial contribution to the low-energy spin Hamiltonian is 
of fourth order, and the operator S is approximated as 

s«E*=iS (n) - 

The unperturbed Hamiltonian, Hq = U\ + U2 + Ub, 
describes localized electrons, and the hopping H tU n acts 
as perturbation. The low-energy subspace of the unper- 
turbed Hamiltonian is spanned by the states in which the 
magnetic ions are singly occupied, and the bridge is dou- 
bly occupied. The lowest-order terms that give rise to a 
nontrivial spin Hamiltonian, in the limit |t|, |P| -C U, are 
of the fourth order in t and P. 

The resulting interaction of the spins includes an 
isotropic exchange of strength J, a Dzyalozhinsky-Moriya 
interaction described by a vector D, and an anisotropic 
exchange term described by a second rank symmetric 
traceless tensor T [601 



H12 — JSi ■ S2 



■D-(Si 



Si-rs 2 . 



(58) 



Quite generally the interaction between two spins up to 
second order in P12 can be represented as an isotropic ex- 
change of rotated spins [HTJ] . However, since the frustra- 
tion in the triangle is strong, it is a good approximation 
to take only the Dzyalozhinsky-Moriya interaction into 
account for the weak spin-orbit coupling, |F*i 2 1 <C |*xa | 
when describing a full molecule. 

In a bond with a single bridge site, the largest possible 
symmetry is G2V We introduce Cartesian coordinates 
with the x-axis pointing from the magnetic center 1 to 
2, y-axis lying in the bond plane and pointing towards 
the bridge site, and the z-axis normal to the bond plane 
(Fig. [SJ. The elements of the C2V symmetry group are 
then rotation R v ^ by 7r about the y-axis, reflection a v in 
the yz plane, and reflection ah in the xy plane. Each of 
these symmetry operations present imposes constraints 
on the parameters of Hb . In the case of localized orbitals 
that remain invariant under the local symmetries of their 
respective sites, the constraints resulting from the R y ^ 



h 
P x ,i 

P z .i 



(57) The a v symmetry implies: 



h 

Px,l 

Py,i 

PzA 



*2, 

— Px,2, 

Py,2> 

-Pz.2- 



Px.2, 



P 



y,2> 



-P, 



and the ah symmetry implies: 

Pi = -P a = pe z 



(59) 
(60) 
(61) 
(62) 



(63) 
(64) 
(65) 
(66) 



(67) 



In the perturbative calculation of the effective spin 
Hamiltonian parameters, these constraints reproduce the 
Dzyalozhinsky-Moriya rules. We do not deal with the 
symmetry of on-site energies t/i,2,b m & n y detail, since 
they do not affect the spin Hamiltonian at this level of 
approximation. 



D. Electric field along y 

In the electric field pointing along the y axis, the point 
group symmetry of the bridge remains C2V , and all of the 
constraints (l60l) - (|67|) hold. The fourth-order Schrieffer- 
Wolff transformation then gives the interaction between 
the spins on magnetic centers of the form (|5"5|) with the 
parameters 



D 



1 



(48i 4 - 40t 2 p 2 z + 3p 



12U 3 
^t Pz (4t 2 -^)e 2 



4), 



(68) 
(69) 



r ' xx — ^yy 2^ ZZ ~ 3C/ 3 ^ (^) 

while all the off-diagonal elements of T vanish. Here, the 
parameters of the Hubbard model satisfy the symmetry 
constraints of the full C2V, and 



h 

Pi 



*2 = t, 

-P2 = p z e z . 



(71) 
(72) 



We have introduced t/ 3 = U c2 (2U c2 - U h2 )(U b i - U h2 + 
{/ C 2) 2 /(4C/ C 2 — Uta), where the on-site repulsions are Ub2 
for the doubly occupied bridge, Ubi for the singly occu- 
pied bridge, and U C 2 for the doubly occupied magnetic 
center. The parameter U describes the energy cost of 
leaving the manifold of states with the minimal energy 
of Coulomb repulsion. We assume that the lowest energy 
charge configuration corresponds to a doubly occupied 
bridge, so that Ub2 < Ubi- 
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In first order, the variations of the spin-Hamiltonian 
parameters resulting from the modification of the Hub- 
bard model parameters, are: 

8 J = ^3 [(48i 3 - 20tp 2 z ) St + (-20t 2 Pz + 3p 3 z ) S Pz ] , 

(73) 

5D Z = A [(\2t 2 Pz - pi) St + (At 3 - M P 2 Z ) 5 Pz ] , (74) 
ST ZZ 16tp z 



ST 



U 3 

sr,. 



w 3 



( Pz St + tS Pz ). (75) 



Electric field modifies the orbitals and therefore the over- 
laps between them, that determine the hopping param- 
eters. We consider the case where the variations St and 
5p z are linear in the field intensity E y : St = KtE y , 
Sp z = Kp„E y . We will not discuss the effect of varia- 
tions in the on-site energies U in any details, since their 
only effect in the fourth order perturbation is a rescaling 
of all the spin Hamiltonian parameters by U 3 / (U + 5U) 3 . 

We stress that these linear modifications of the hop- 
ping parameters are characteristic for the C2v symmetry. 
If the electric field is oriented differently and thus lowers 
the system symmetry (see below) first-order increments 
are not allowed, and the spin-electric coupling is at least 
a second order effect in the electric field. The modifica- 
tion of the orbitals includes the energy scale of splitting 
of the atomic orbitals in the molecular field. We have 
assumed earlier that the splitting of the orbitals local- 
ized on the magnetic centers is large, and the dominant 
source of the spin-electric coupling is the modification of 
the bridge orbital. Therefore, the key criterion for strong 
spin-electric coupling is the presence of bridge orbitals 
that are weakly split in the molecular field. If, in addi- 
tion, we assume that the modification is a property of 
the bond alone, and not of the entire molecule, the k pa- 
rameters can be determined in an ab-initio calculations 
on a smaller collection of atoms. 

In the limit of weak spin-orbit coupling, \t\ 3> \p z \, the 
main effect of the electric fields is a change of J, leading 
to our symmetry-based results, see EqUHl In particular, 
the d parameter of the symmetry analysis is: 



d= — [(48i 3 - 20tp 2 z ) n t + 



-2tit 2 p z + 3p 3 z ) k pz ] 



(76) 



In this case, the Dzyalozhinsky-Moriya vector D is 
constraint to point in the z direction, D = De z . The 
model suggests that the dominant effect of the electric 



field in the molecules with dominant Heisenberg exchange 
( J ^> |D|) is modification of the isotropic exchange con- 
stants J, and 



\SD[ 

\SJ\ 



\J\ 



(77) 



so that the modification of the Dzyalozhinsky-Moriya 
vector D — ► D + is weaker. However, in the molecules 
in which the modifications of J are inefficient in inducing 
the spin-electric coupling, as for example in the spin- 1/2 
pentagon, the modifications of D may eventually provide 
the main contribution to the spin-electric coupling. 

Electric field pointing in a generic direction breaks the 
C2v symmetry of the bridge, and allows further modifica- 
tion of the Hubbard and spin Hamiltonian parameters, 
that do not obey all the symmetry constraints in Eqs. 
(|59"]l - (joTf . With the relaxed constraints, both the direc- 
tion and intensity of Pi, 2, as well as the spin-independent 
hoppings ti_2 become field-dependent. This observation 
can be used in the search for molecules that show strong 
spin-electric coupling. The energy cost of changing the 
distance between the localized orbitals may be signifi- 
cantly higher than the cost of modifying the shape of the 
bridge orbital. In order to investigate this dependence, 
we study the effective spin Hamiltonian description of a 
bridge with all possible residual symmetries. 



1. Residual a v symmetry 

An electric field E = Ee z normal to the bond's plane 
reduces the initial C21, symmetry down to {E,a v }. This 
reduction of the symmetry also happens when a molecule 
is deposited on the surface parallel to the bond plane. 
While the constraints in Eq. [57] hold, this reduction 
of symmetry implies the appearance of nonzero in-plane 
components of Pi, 2- We parameterize the most general 
Hubbard model parameters £1,2, Pi. 2 consistent with the 
symmetry as 



t\= t2 = t, 

Pl,x = P2,x = Pxy COS (f>, 
P\,y = ~P2,y = Pxy Sine 
Pl,z = -P2,z = Pz- 



(78) 
(79) 
(80) 
(81) 



The effective low energy spin Hamiltonian, derived by 
Schrieffer- Wolff transformation up to fourth order in t/U , 
and |P|/[/ is given by ([55]) . with the non-zero parameters 
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J = 



Dy 

D z 
r 

x XX 

V 

1 yy 
v 

yz 



[pt y - 2p 2 xyP 2 z + 3pt - St 2 (p 2 xy + 5p 2 z ) + 48< 4 - 8p 2 xy (p\ - U 2 ) cos 20 + 2p% cos < 
-jjj(Pz cos + 2i sin 0) (-p 2 + 4i 2 + p 2 xy cos 20) , 
- 7^p (4tp z - pl y sin 20) (p 2 - 4i 2 - p 2 xy cos 20) , 

: -Qlp [Ply (! " cos2 ^) + 2 P'] t 8 * 2 + C 1 + cos2 ^)] > 



12C/ 3 



{ -Pty + SplyPl + 32t 2 (pl y ~pl)+ pI v [8 (p 2 z - it 2 ) cos 20 + p 2 X y cos 40 + 48tp z sin 20] } , 



= T zy = ^ (p a cos + 2* sin 0) (-Atp z + p 2 xy sin 2c, 



■y 

= -T 



(82) 
(83) 
(84) 
(85) 
(86) 
(87) 



- yy 



In the lowest order in spin-orbit coupling the spin in- 
teraction consists of the isotropic exchange with J rj 
4t 4 /U 3 , and the Dzyalozhinsky-Moriya interaction with 
D w -8^^^ sin0e a + p z e z )/[/ 3 . 



As a matter of principle, the spin-orbit coupling medi- 
ated hopping P does not have to be much weaker than the 
spin-independent hopping t. In this case, all the nonzero 
terms in Eqs. (|82|) — (|88| are of comparable size, and the 
variation of spin Hamiltonian with the angle becomes 
significant. Note that the angle describes the directions 
of spin-orbit coupling induced hopping parameters Pi. 2, 
and that it is not directly connected to the bond angle 
between the magnetic sites and the bridge site. How- 
ever, for the bridge orbital without azimuthal symmetry, 
the angle does depend on the bond angle. For the 
molecules in which the full symmetry allows only for the 
spin-electric coupling mediated by the spin-orbit interac- 
tion, this effect is important. 



suiting change in the spin-Hamiltonian parameters reads: 

+ Pz0 (-20t 2 + 3 P 2 zQ ) 6 Pz ] , 

SD y = ~jp ( 4 <o -Pto) (2fosin0+p zO cos0) 6p xy , 

(90) 

6D, = [p zQ (l2t 2 -p 2 z0 )5t (91) 

+t (4i§ - 3 P 2 z0 ) S Pz ] , 
ST XX = 5Tyy = --ST ZZ = (92) 

- T^toPzO (PzoSt + t Sp z ) , 

4 

ST yz = ST zy = -— ^t p z0 (2t sin + p z o cos (f>)5p xy . 

(93) 



With these assumptions, the dependence of the effec- 
tive spin Hamiltonian on p xy suggests that the strength 
of induced in-plane Dzyalozhinsky-Moriya vector will be 
sensitive to the angle that is determined by the angular 
dependence of the bridge- and magnetic center orbitals. 
In turn, for a fixed symmetry of the bridge orbital, this 
dependence directly translates into the dependence of the 
spin-electric coupling constant on the bridge bond angle. 



In the presence of electric field E = Ee z , the hopping 
parameters will change from their initial values, that sat- 
isfy the constraints implied by the C2 V symmetry, into a 
set of values that satisfy those implied by a v only. The re- 



The cr^-symmetric variations of Hubbard parameters 
occur when an external electric field is applied along 
the z direction to a C2V symmetric bond. Again, the 
variations of the parameters is generically linear in the 
field strength, St = Kt,<n>Ez, $Px y = K P x y ,crvE Zl Sp z = 
i^pz,crvE z , where the k parameters depend on the mod- 
ification of the bridge orbital in the electric field. As 
opposed to the case of the field along y direction that 
maintains the bonds C2V symmetry, the k parameters 
for the field along z axis vanish in zero field, since the 
z-component of a vector has no matrix elements between 
the relevant C2„-symmetric states. The linear expansion 
is valid when the field is strong enough to distort the 
bridge orbital. Alternatively, the expansion is valid for 
a bond with lower symmetry in zero electric field, e. g., 
when the bond is close to a surface. 
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2. Residual ah symmetry 

In an electric field that lies in plane of the bond, with 
E || x), the only residual symmetry transformation is the 
reflection about the xy plane (ah). Within this reduced 
symmetry, the two magnetic sites are no longer equiv- 
alent, but the spin-dependent hopping parameters Pi. 2 
still point along the z axis: 



ii ^ t 2 , Pi = Pie z ^ p 2 e z = P 2 



(94) 



In the fourth order in hopping t, P, the resulting low 
energy spin Hamiltonian is again given by Eq. I58[ with 
the following non-zero coupling constants: 



J 



1 



[32t 1 t 2 Pl z P2z 



12U 3 

-4 (t\ v \ z + t\ v \ z ) + m\t\ + 3 P LpL] , 

D = -^3 (hp2z - hPiz) (4iii 2 + PizP2z) e z , 

\x = T TO = ^- = -7^p [hP2z ~ t2Plz f , 



(95) 

(96) 
(97) 



Similarly to the case of full C21, symmetry, the 
spin Hamiltonian consists of the isotopic exchange J, 
Dzyalozhinsky-Moriya vector D = D z z normal to the 
bond plane, and diagonal tensor T isotropic in the bond 
plane (T xx = F yy ). We stress that the dependence of the 
effective spin Hamiltonian parameters on those entering 
the spin Hubbard Hamiltonian is different for these two 
symmetries, and so is the response to the applied elec- 
tric field. On one hand, the C2V preserving electric field 
induces the transitions in the lowest energy multiplet in 
the lowest order. On the other hand, the electric field 
that reduces the bond symmetry to {E,ah} does not al- 
ter the coupling of spins in the lowest order, since the 
deformation of the molecule requires some coupling to 
the field. 

As in previous case, we expand the ah symmetric spin 
Hamiltonian around the C2V symmetric case. We intro- 
duce a perturbation of the parameters Hubbard parame- 
ters in the electric field consistent with the residual sym- 
metry: tx = t + Sti, t 2 = t + St 2l piz = Pzo + Spiz, 
P2 Z = —Pzo + Sp2 Z - As a consequence, the spin Hamilto- 
nian parameters are incremented by: 



6J = 



1 

61P 



[4t (I2tg - 5^o) + St 2 ) 



+ PzQ (-20t§ + 3p 2 z0 ) (S Plz - S P2z )] , 
SD Z =— [ Pz0 t a (12t Q ~ Pzo) (oh + 6t 2 ) 
+t (4tg - 3p2 Q ) (fipu - S P2z )] , 

sr zz 



5T XX =5T 



8 
3LP 



hPzO [PzO (Sti +St 2 ) + t (dpi 



(99) 

(100) 
Sp 2z )} ■ 



As for the case of a v residual symmetry, there is no spin- 
electric effect of the first order in electric field, and the 
crucial condition for coupling to the electric field in this 
direction is weak splitting of the bridge orbitals in the 
molecular field. 



Residual R y .-w symmetry 



Reduction of the symmetry of the bond, from the full 
C2„ to the group {E, Ry,-*}, does not occur for any vector 
perturbation. In terms of electric fields, this reduction 
of the symmetry would correspond to an inhomogeneous 
electric field that points in the e z direction at the position 
of one of the magnetic centers, and in the — e z direction 
at the position of the other. This symmetry breaking can 
also happen when the localized orbitals on the magnetic 
centers have lobes of opposite signs extending in the z- 
direction, and oriented opposite to each other. 



The most general Hubbard model parameters consis- 
tent with the residual symmetry are 



t\ = *2 = t, 
Plx = -P2x = Pxy COS ( 
Ply = P2y = Pxy SU1 (f), 
Plz = P 2z = P Z - 



(101) 
(102) 
(103) 
(104) 



After the fourth-order Schrieffer- Wolff transformation, 
the effective low-energy spin Hamiltonian has the form 
(|58[1 with nonzero parameters 
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J = 



1 



{ Pxy - 2p 2 xyP 2 z + 3p 4 z - 8i 2 {pl y + 5p 2 z ) + 48* 2 + Sp 2 xy {pi - At 2 ) cos 20 + 2 Pxy cos 40) , 



12C/ 3 

= jpPxy (-2* cos + p z sin 0) (p 2 - 4* 2 + p 2 xy cos 20) , 
1 

1 



D z = -^fj3 {^Vz + Ply sin 20) (p 2 - 4t 2 + p 2 xy cos 2 
r = 

12 LP x 

1 



" " 12Z73 (~ p * v + 8p2xyp2z + m2 ( p * y ~ p ^ + p2y ^ 8 ~ 4 ^ cos 2(/) + p2y cos 40 ~ 48 ^ z sin 2 ^ ) 



r — r — 



xz = jj^Pxy (2* cos 4>~p z sin 0) (4*p z + p 2 xy sin 20) , 



T yy = ^773 (*4 (! + cos2 ^) + 2 p 2 *) (pIv (- 1 + cos2( « ~ 8 ^) > 



6U 3 

r„ = -r 



(105) 
(106) 
(107) 
(108) 
(109) 
(110) 



r ^ = -7TT73 + 2 PlyPl + 8<2 (*4 " A pD 



-p 2 (-2(p 2 -4* 2 )cos- J 



20 + p 2 a cos 40 - 24tp z sin 20) ) 

(111) 



r 



The expansion from the C2,, symmetric case gives (see 
the discussion of the a v residual symmetry in Subsection 
11111)11: 

5 J =^3 [4to (12* 2 , - 5 P 2 z0 ) St + Pz0 {-20t 2 + 3p 2 ) 5 Pz ] 



(112) 



5D X = — (4*o - p 2 ) (2*o cos 0o - p z0 sin O ) <5pxj/, 



(113) 



<5D* =— [p z0 (I2£g - p 2 ) St + t {At 2 - 3p 2 ) Sp z0 ] 



ST XX =STyy = -^Y z 



16 
3LP~ 



P 2 o*o (p 2 o^ + t Sp z ) , 



(115) 



<5r 2X =ST XZ = jj^t p zQ (2*0 cos 0o - p z0 sin0 o ) <5p X2/ . 

(116) 

As in the case of a v symmetry, the resulting interac- 
tion of the spins on magnetic centers becomes dependent 
on the angle between the two P parameters. This de- 
pendence is pronounced in the case of strong spin-orbit 
coupling and can lead to the dependence of spin-electric 
effects on both the geometry of the bond and the shape 
of the bridge orbital. 



E. Bond modification and symmetries 

Spin-electric coupling induced by the superexchange 
through bridge atoms depends on the symmetry of the 
bridge and the direction of the electric field. This symme- 
try reflects on the resulting coupling of spins in an MN. In 
this subsection, we combine the results of the Hubbard 
model study of the individual bonds with the previous 
symmetry considerations, and provide rough estimates 



of the most promising spin-electric coupling mechanism 
in the triangular and pentagonal molecules. 

The spin-electric coupling via superexchange is most 
sensitive to the electric fields that does not break the 
initial C2„ local symmetry of the bond. This symmetry 
corresponds to the electric field that lies in the plane 
of the molecule and normal to the bond. All the other 
couplings require modification of the bridge orbitals, and 
are suppressed by a factor d\E\/Ud, where U& is on-site 
repulsion on the bridge. Assuming that this repulsion is 
strong, we can model the spin electric coupling as a set of 
modifications of the spin interactions SHjj + i between the 



(114) neighboring magnetic centers, with oc \E 



bond I 



where E^ nd is the projection of the electric field normal 
to the bond and lying in the molecule's plane. 

In the triangle, the strongest effects of electric 
field is modification of exchange couplings 5Jjj+± = 
5 Jo cos (2j7r/3 + #0), where the angle 80 describes the ori- 
entation of the in-plane component of the electric field, 
and 5 Jo is a molecule-dependent constant. This modifi- 
cation leads to a specific coupling of the in-plane compo- 



nents of chirality to the electric field H' t 



off 



rfE'-Cn, see 



([7]) . Other types of coupling are suppressed either due to 
weaker influence of electric field on the bonds, or due to 
the symmetry of the molecule. If the spin-electric cou- 
pling is mediated by the spin-orbit interaction, the sup- 
pression is by a factor of the order |D|/J, and if the cou- 
pling is mediated by electric field, the suppression factor 
is d\E\/J. Assuming the simplest case, the modification 
of exchange coupling is the most promising mechanism 
for spin-electric coupling in triangular molecules. 

In the pentagons, the modification of spin-spin in- 
teraction SHjj + i preferred by the superexchange mech- 
anism is inefficient in inducing the spin-electric cou- 
pling of the molecule. The pattern 5Jjj+i of exchange 
coupling constants induced by an external electric field 
does not couple the states within the lowest energy 
manifold. In order to couple the spins in the pen- 
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tagon to an external field, another mechanism is needed. 
The modification of the Dzyalozhinsky-Moriya vectors 
5Dj z — SD z o cos (2j7r/5 + 8q), where 5D z o is a molecule- 
dependent constant, and 8q describes the orientation of 
the in-plane component of the electric field, is preferred 
by the superexchange bridge model. In the symmetry 
analysis, we have found that this form of modification 
of spin-orbit coupling does not induce spin-electric cou- 
pling. The same applies to the modifications of in-plane 
components T)j xy . The main effect that gives rise to 
spin-electric coupling is the modification of the exchange 
interactions 8Jjj+\ in the presence of the original spin- 
orbit interaction Djj+\ tZ . Compared to a triangle com- 
posed out of identical bonds, this interaction will be 
weaker by a factor of |Djj+i|/Jjj+i. 

In summary, within our model of the superexchange- 
mediated spin-electric coupling, the most promising can- 
didates for the spin manipulation via electric field are tri- 
angular molecules. In pentagons, the best candidates are 
molecules with strong spin-orbit interaction, and weakly 
split bridge orbitals. 

IV. EXPERIMENTAL SIGNATURES OF THE 
SPIN-ELECTRIC COUPLING 

Coherent quantum control of spins in an MN using 
electric fields can be achieved by resonant driving of 
the transitions between the chirality eigenstates [4l|. At 
present, however, little is known about the effects of elec- 
tric fields on the spin states of molecular magnets. As a 
preliminary step, it is useful to identify possible signa- 
tures of such a coupling that are observable in the exper- 
iments routinely used to characterize these systems. 

In this section, we study the ways in which the spin- 
electric coupling can be detected in electron spin reso- 
nance (ESR), in nuclear magnetic resonance (NMR), and 
in the thermodynamic measurement of an MN. 

A. Electron spin resonance 

Electron spin resonance (ESR) investigates transitions 
between states belonging to a given S multiplet and hav- 
ing different spin projections M along the magnetic field 
direction [6lj ]. This technique provides information on 
the anisotropies of the spin system, as well as on the 
chemical environment, and the spin dynamics [29| . In the 
following, we show how the effects of an external electric 
field can show up in the ESR spectra of antiferromag- 
netic spin rings by affecting both the frequency and the 
oscillator strength of the transitions. 

1. Triangle of s = 1/2 spins 

We start by considering the simplest case of interest, 
namely that of a triangle of s = 1/2 spins with D^h sym- 



metry. The lowest energy eigenstates of the spin triangle, 
given in Eq. © form an S = 1/2. The effective Hamilto- 
nian H e g of the molecule in the presence of electric and 
magnetic field, and acting within this quadruplet is given 
by Eq. (15]). 

We first consider the case of a static magnetic field per- 
pendicular to the molecule's plane (B || z). The eigen- 
values of H e g are then given by: 

K = <y[B + a(Al +£Y /2 l (H7) 

where £ = d\B x z|, B = HB^gjBl + g\ B\ , a = ±1/2 

is the eigenvalue of S z , a = ±1 is the the eigenstate 
chirality in the limit of vanishing electric field: |A^)f = o = 
(— a) cr ~ 1 / 2 |o;, a). In the presence of electric field, the 
eigenstates read: 

|A«) = {2a[A so + a(£ 2 + A 2 o y/ 2 }\ + l,a) 

+ £e- l8 \ - l,a)}/D a , (118) 

where D a = {£ 2 + [A so + a(£ 2 + A 2 ^) 1 / 2 } 2 } 1 / 2 . 

Electron spin resonance induces transitions between 
such eigenstates. The transition amplitudes are given by 
the absolute values of matrix elements of x-component of 
the total spin, taken between the states that the transi- 
tion connects, 

(A" 1/2 |S,|A;? /2 ) = -S'/D+'D- 1 (119) 

, x a \q \ \ a \ - AsotAso + a^ + Ajy 1 / 2 ] 
(A_ 1/2 |/b x |A +1/2 ; - T--T2 (-t2U) 

The corresponding frequencies are given by: 

A +i/ 2 - Al? /2 = B (121) 
K1/2 ~ A-i/ 2 = B + a{£ 2 + A 2 so ) Va. (i 22 ) 

As an illustrative example, we plot the frequencies and 
amplitudes of the ESR transitions as a function of the 
electric field (Fig. [6]). While for £ = 0, these transitions 
can only take place between states of equal C z (red and 
green symbols online, transitions with the larger ampli- 
tude at low fields, in the figure and in the inset), the elec- 
tric field mixes states of opposite chirality, thus transfer- 
ring oscillator strength to two further transitions, whose 
frequencies are independent of £ (blue symbols online, 
constant frequency transition in the figure). In the limit 
dE S> D z , the eigenstates of the spin Hamiltonian tend 
to coincide with those of Sf 2 , and ESR transitions take 
place between states of equal Si2- While the eigenstates 
depend on the in-plane orientation of the electric field, 
no such dependence is present in the frequencies and os- 
cillator strength of the ESR transitions. Besides, these 
quantities are independent of the exchange coupling J, 
and depend on the value of the applied magnetic field 
only through an additive constant (cjq). 

The dependence of the ESR spectrum on the applied 
electric field is qualitatively different if the static mag- 
netic field is applied in-plane (e.g., B || x and the oscil- 
lating field oriented along z). In this case, the eigenvalues 
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FIG. 6: (color online) Energy (cj) of the ESR transitions in a 
triangle of s = 1/2 spins as a function of the applied electric 
field E that lies in the molecule's plane, so that d|E| = dE = 
£ . The magnetic field is B || z and ujq — gfisB, see Eqs. 
(|121|) and (| 122f) . The diameter of the circles is proportional 
to the transition amplitudes | (a\S x \ot') |, Eqs. (|119|l and (|120|1 • 
Here, |a) are the eigenstates of H in the lowest energy 5=1/2 
multiplet. Inset: Eigenvalues (in units of D z ) as a function of 
£ = d|E|||, in units 3D z /4. 



FIG. 7: Energy (u) of the ESR transitions in a triangle of 
s = 1/2 spins as a function of the applied in-plane electric 
field E, so that d[E\ = dE = £ , and in the presence of the 
in-plane magnetic field B || x. The diameter of the circles is 
proportional to \(a\S z \a')\, Eqs. (fmj) and {122]). The states 
| a) are the eigenstates of H in the lowest S — 1/2 multiplet. 
Inset: Eigenvalues (in units of D z ) as a function of d|Eii | = £, 
in units 3D z /4. 



of -f/eff are: 

< = aa[A 2 + (£ + aB) 2 }^ 2 , (123) 

whore (7 = ±1/2 is the value of (S x ) in the limit of large 
magnetic field (B ^ £,Aso) and a = ±1. The corre- 
sponding eigenstates read: 

K) = {e w (Aso + ^)[| + l,+l/2)-|-l-l/2}] 

+ (B + a£)[\ + l -1/2) - |-l,+l/2)]}/^124) 

where 

D a a = V2[(A so +^) 2 + (B + a£) 2 }^ 2 . (125) 

The expectation values of the total spin along the mag- 
netic field for each of the above eigenstates are given by 
the following expressions 

{tf\S x \tf) = 2[(A so + ^)(£-a£)]/ra 2 ,(126) 

which are independent of the in-plane direction of the 
electric field. The ESR transitions between such eigen- 
states induced by a magnetic field that oscillates along 
the z direction are given by the expressions: 

, a]q | -m _ (A S o + /ig)(A so + M -") + (g 2 -g 2 ) 

(f£\SM) = 0. ^ " (127) 

Therefore, the application of the electric field shifts the 
energy of the transitions between states of opposite a, 



thus removing their degeneracy; however, unlike the case 
B || z, it does not increase the number of allowed transi- 
tions. 

In the case of tilted magnetic fields, the dependence of 
the ESR spectrum on the applied electric field presents 
qualitatively different features (Fig. [?])• In particular, 
the spectrum is dominated by two pairs of degenerate 
transitions that anticross as a function of the electric 
field. Away from the anticrossing, the transitions with 
the largest oscillator strength display frequency depen- 
dence on the electric field. 



2. Pentagons of s = 1/2 spins 

Triangles of s — 3/2 spins (not shown here) display the 
same qualitative behavior as the one discussed above. In 
contrast, chains including an odd number N > 3 spins 
behave differently. This is mainly due to the fact that 
the spin-electric coupling 8H does not couple directly 
the four eigenstates of H belonging to the lowest 5 = 1/2 
multiplet: such coupling only takes place through mix- 
ing with the higher S = 1/2 multiplet. As a consequence, 
the effects of the spin-electric coupling tend to be weaker 
as compared to the case of the triangle; besides, unlike 
the above case of the spin triangle, they depend on the 
exchange coupling J. Illustrative numerical results are 
shown in Figs. [S] and H] for the cases of a perpendicu- 
lar and in-plane magnetic field, respectively. In the for- 
mer case, both the frequencies and amplitude of the ESR 
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FIG. 8: Energy (w) of the ESR transitions in a pentagon of 
s — 1/2 spins as a function of the electric field applied in 
the molecule's plane d|E| — dE — £. The Zeeman splitting, 
u)o — gi-iBB is set by the magnetic field B || z, orthogo- 
nal to the molecule's plane. The considered transitions are 
those between eigenstates (\a)) belonging to the 3 = 1/2 
multiplet of the spin Hamiltonian (figure inset). Unlike the 
case of the spin triangle, these are coupled to each other by 
the electric field via eigenstates belonging to other multiplets, 
and therefore depends also on the exchange constant J (here 
J/Aso = 100). The diameter of the circles is proportional to 
| {a\S x \<x') |, and therefore to the transition amplitude. 
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FIG. 9: Energy (lu) of the ESR transitions in a pentagon of 
s — 1/2 spins as a function of the applied in-plane electric 
field E, so that d|E| = dE — £. The Zeeman splitting is set 
by an in-plane magnetic field B || x, and ljo = gfisB. The 
considered transitions are those between eigenstates (|a)) be- 
longing to the S = 1/2 multiplet of the spin Hamiltonian 
(figure inset). Unlike the case of the spin triangle, these are 
coupled to each other by the electric field via eigenstates be- 
longing to other multiplets, and therefore depends also on the 
exchange constant J (here J/Aso = 100). The diameter of 
the circles is proportional to Ka^S^c/}!, and therefore to the 
transition amplitude. 



transitions are hardly affected by the electric field, in the 
same range of physical parameters considered in Fig. 
In the case of an in-plane magnetic field, instead, a rela- 
tively small shift in the transition energies is accompanied 
by a strong transfer of the oscillator strength, for values of 
the spin-electric coupling exceeding the Dzyalozhinsky- 
Moriya coupling constant. 



B. Nuclear magnetic resonance 

The spin-electric Hamiltonian 5Hq modifies non uni- 
formly the super-exchange couplings between neighbor- 
ing spins. This might not affect the projection of the 
total spin (as in the case B || z, see above), but it gen- 
erally affects the moment distribution within the spin 
chain. Such effect can be investigated through exper- 
imental techniques that act as local probes in molec- 
ular nanomagnets, such as nuclear magnetic resonance 
(NMR) [62| or x-ray absorption (§1}. In NMR, the ex- 
pectation value of a given spin within the cluster can be 
inferred through the frequency shift induced on the tran- 
sitions of the corresponding nucleus. The shift in the 
nuclear resonance frequency for the nucleus of the i-th 
magnetic ion is 1S.V = r )A(s z ^), where A is the contact 
hyperfine interaction constant at the nuclear site. The 
constant of proportionality A depends on the spin den- 



sity at the position of the nucleus, and can be extracted 
from the experiment by considering the polarized ground 
state M = S at high magnetic fields [62[ . As in the case of 
ESR, the dependence of the NMR spectra on the applied 
electric field qualitatively depends on the orientation of 
the static magnetic field B with respect to the molecule. 
Unlike the case of ESR, however, it also depends on the 
in-plane orientation of the electric field, i.e. on the way 
in which the E breaks the symmetry of the molecule. 



1. Spin triangles 

Let us start by considering a spin s = 1/2 triangle, with 
a magnetic field applied perpendicular to the molecule 
plane (B || z). In this case, the distribution of the spin 
projection along z is given by the following expression: 

(A?|«i,,|A°)=<7/3 + /»(f)cos[fl + 7r(5/3-i)], (128) 
where 

^[Aso + cKAlo+f 2 ) 172 ] 



(129) 



Here, the expressions of the eigenstates |A") and of D a 
are given in Subsection IIV Al For E = 0, the three 
spins are equivalent and (A" x /2 = ±1/6. 
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FIG. 10: (Color online) Expectation values of the z- 
component of s — 1/2 spins in a triangular molecule as a 
function of applied electric field. The magnetic field is per- 
pendicular to the ring plane (B || z); the electric field is par- 
allel and perpendicular to ri2 in the upper and lower panel, 
respectively. In the electric field along one of the bonds (lower 
panel), the spins that lie on that bond have the same out-of- 
plane projections. The shadings (colors online) denote the 
different spins. 



If the electric field is finite and oriented along one of the 
triangle sides (e.g., E || ri 2 , corresponding to 6 = 0), 
then expectation values along z of spins 1 and 2 un- 
dergo opposite shifts, whereas that of spin 3 is left un- 
changed: A E (s M ) = -A E (s 2 ,z), where A E (s i;2 ) = 
( s j,z)e — (si,z)e=o- This is shown in Fig. [10] for the 
ground state of the spin Hamiltonian, but the above re- 
lations hold for any of the four eigenstates |A") belong- 
ing to the S = 1/2 quadruplet. If the NMR frequency 
shifts Afi are larger than the corresponding line widths, 
the single line at E — splits into three equispaced 
lines, with intensity ratios 1:1:1. If, instead, the elec- 
tric field is applied along a symmetry plane of the tri- 
angle (e.g., E _L ri2, corresponding to 9 = ir/2), spins 
1 and 2 remain equivalent and their magnetic moments 
display the same electric field dependence, while the shift 
of the third one is opposite in sign and twice as large in 
absolute value: A E (si :Z ) = A E (s 2 , z ) = -A E (s 3 , z )/3. 
The intensity ratios of the two NMR lines are, corre- 
spondingly, 1:2. The expectation values for the remain- 
ing eigenstates can be derived by the following equa- 



/2l s i^|-^-l/2/ 



-A E (A« 1/2 | S 



|A Q 



1/2/ 



tions: A E (A^ 1 

and (Xl\8 i>s \Xl) = -A E (A- 1 |s l , 2 |A- 1 ). Therefore, at fi- 
nite temperature, the shifts in the expectation values of 
the three spins are given by: 



A E (s i; 



(130) 



A E (A^ /2 K z |A+i /2 ) J2 a cosh(j^) 
If the field is oriented along the molecule plane (B || x), 



FIG. 11: (Color online) Expectation values of the z- 
component of s — 3/2 spins in a triangular molecule as a 
function of applied electric field. The magnetic field is per- 
pendicular to the ring plane (B || z); the electric field is par- 
allel and perpendicular to ri2 in the upper and lower panel, 
respectively. The shadings (colors online) denote the different 
spins. 



the expectation value of the three spins corresponding to 
each of the eigenstates are given by the following expres- 
sions: 

- + (1/3) cos(0 - 2*7t/3), (131) 



where 



2 (Ago 
3 



fj£)(B + c£) 



(132) 



If the magnetic field is parallel to the triangle plane, 
the in-plane electric field can modify the total spin ex- 
pectation value along B. The changes that E induces 
in the magnetization distribution within the triangle at 
zero temperature are less varied than in the previous case 
(Fig. [TU]). In fact, the magnitude of the A E (si )Z ) is much 
smaller, and all the spins undergo shifts of equal sign and 
slope. The NMR line, which is slitted into three lines al- 
ready for E = 0, is rigidly by the applied electric field. 

If the triangle is formed by half-integer spins s > 1/2, 
an analogous dependence of the expectation values (sj, z ) 
on the electric field is found. As an illustrative example, 
we report in Fig. [11] the case of s = 3/2. 



2. Pentagon of s = 1/2 spins 

Spin chains consisting of an odd number of half- integer 
spins present analogous behaviors, but also meaningful 
differences with respect to the case of the spin triangle. 
In particular, the spin-electric Hamiltonian SHq does not 
couple states belonging to the lowest S = 1/2 quadruplet 
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energy scale in the spin Hamiltonian, and being mainly 
interested in the low-temperature limit, we restrict our- 
selves to the 5 = 1/2 quadruplet, and use the effective 
Hamiltonian H e g in Eq. (|2"Tj) . 

The eigenenergies of the lowest 5 = 1/2 sector in the 
presence of electric and magnetic fields are 



with B = HB\Jg 
(i3£) 2 )] 1/4 , and B z 



B 2 



A 2 



£ 2 + 2 1 E 2 1 (133) 
+ g\Hl, E = l(B z A so ) 2 + 

Hb9\\H\\. Note that these ener- 



gies are the generalization of the ones in the previous 
section, which were valid for in-plane magnetic field only. 
The partition function for N identical and non-interacting 
molecules is Z = Zy , with Z\ = J2 a 7 ex P ( — flEa,j) 
being the partition function for one molecule, and (3 — 
l/(k-QT). The free energy reads 



F = -1/pinZ = -Nk B Tln 



2^cosh(/3.E 7 ) 



, (134) 



FIG. 12: (Color online) Expectation values of the z- 
component of s = 1/2 spins in a pentagon as a function of 
applied electric field. The magnetic field is perpendicular to 
the ring plane (B || z); the electric field is parallel (6 = 0) 
and perpendicular (6 = 7r/2) to ri2 in the upper and lower 
panel, respectively. The shadings (colors online) denote the 
different spins. 

directly (i.e., matrix elements (i\SHo\j) = for i,j < 
4); these couplings are mediated by states belonging to 
higher 5 = 1/2 multiplets, that are higher in energy by 
a quantity ~ J. Therefore, the effect of the electric field 
tends to be significantly smaller than in the case of a 
triangle with equal D z and S (see Fig. [T2l , and depends 
also on the exchange coupling J. 

C. Magnetization, Polarization, and 
Susceptibilities 

The spin-electric coupling shifts the energy eigenvalues 
of the nanomagnet, thus affecting thermodynamic quan- 
tities, such as magnetization, polarization and suscepti- 
bilities. In the following, we compute these quantities in 
the case of the s = 1/2 spin triangle as a function of the 
applied magnetic and electric fields. Under the realistic 
assumption that the exchange splitting J is the largest 



with E n = E 1 / 2tl - From this, we can derive different 
thermodynamic quantities like the magnetization Mi = 
—dF/dHi, the electric polarization Pj = —dF/dEi, 
the heat capacity C = —d/dT(dln(Z)/d/3), and the 
corresponding susceptibilities: XEiE, = dPi/dEj = 
d 2 F/dEidEj - the electric susceptibility, XHiHj = 
dMi/dHj = d 2 F/dHidHj - the spin susceptibility, and 
XEiHj = dPi/dMj = d 2 F/dEidHj - the spin-electric 
susceptibility. For the electric polarization components 
Pi we get 

_ Nd£i -sp sinh(/3i; 7 ) 

while for the magnetization components Mj we get 

= NgiUBBi sinh (J3&y) 

1 2E 7=±1 cosh(/?£ 7 ) 7 ^ i E 1 

x (l + 7 ^±£!) l (136) 

where again i — x, y. Making use of the above expres- 
sions, we can obtain the above defined susceptibilities 
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FIG. 13: Electric polarization P x (x component) in Eq. (|135p 
as a function of the magnetic field in x direction. The three 
lines correspond to various values of an additional external 
electric field in the z direction. The plot is for the temperature 
k B T = O.OOlAso, and the electric field dE x = 0.1A S o- 
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FIG. 14: Electric polarization P x (x component) in Eq. (|135p 
as a function of the magnetic field in z direction. The three 
lines correspond to various values of the external magnetic 
field in the x direction. The plot is for the temperature feT = 
O.OOlAso, and the electric field dE x = 0.1A S o- 



The polarization P, magnetization M, and susceptibili- 
ties x, Eq- (|135l) - Eq. (|139l) . all depend on the spin- 
electric coupling constant d. In the following, we analyze 
the details of this dependence and identify the conditions 
suitable for extracting the value of d from the measurable 
quantities. 



1. Polarization and magnetization 

The in-plane polarization of the molecule as a function 
of the magnetic field is illustrated in Fig. [13] and Fig. 
1141 The polarization is a growing function of the mag- 
netic field strength, and it gets reduced by the normal 
component of the field. 

The low-temperature, k^T <C Ago, thermodynamic 
properties of a molecule with spin-electric coupling show 



a simple dependence on the strength of external electric 
and magnetic fields in the special cases of in-plane and 
out-of plane magnetic field. We focus only on effects in 
leading orders in electric field under the realistic assump- 
tion that the electric dipole splitting is small compared 
to the SO splitting, i.e. £ <C Ago- Also, we analyze 
two limiting cases: (i) k^T <C £ , i.e. low-temperature 
regime, and (ii) fcsT 3> £ , i.e. high temperature regime. 
However, we assume all temperatures (in both regimes) 
to satisfy ksT -C Ago so that the spin-orbit split levels 
are well resolved. In the first case (i), we obtain for the 
polarization 



Pi 



nd£, L B 

ndA§ £i 
4A| 



for £ <C 2? 



for £ > B , 



(140) 
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while for the second situation (ii) we obtain 



Pi 



ndAg £i 



B 2 

A 2 
^so 



(3A t 



(141) 



with A B = ^B 2 + A| Q and n = N/V the density of 
molecules in the crystal. We see that, for low tempera- 
tures, the electric polarization Pi ranges from being in- 
dependent of the magnitude of the electric field « 6), 
to a linear dependence on the applied electric field E for 
large fields (£ >• B). Also, the polarization is strongly 
dependent on the magnetic field (linear in B) for low 
i?-ficlds, thus implying strong magneto-electric response. 

We now switch to the other special case, namely when 
the external magnetic field is applied perpendicularly to 
the spin triangles. The electric polarization now reads 



P 



nd£i 
4A7 



tanh (pAs), 



(142) 
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FIG. 15: In-plane magnetization M x in ^-direction in Eq. 
(|136[) as a function of the electric field E x in ^-direction. The 
three lines correspond to a fixed value of an additional mag- 
netic field in the z-direction. The assumed temperature is 
ksT = O.OOlAso, while in the inset it is at higher tempera- 
ture k B T = O.lAso- 



with As = \J Ag + £ 2 . The polarization Pi does not 
depend on the magnetic field B, and there are no spin- 
electric effects present for this particular case. 

Our results suggest that the spin-electric coupling can 
be detected by measuring the polarization of the crys- 
tal of triangular single molecule antiferromagnets that 
lie in parallel planes in the in-plane electric and magnetic 
fields. 

The out-of plane component M z of the molecule's mag- 
netization is rather insensitive to the electric fields, since 
any effect of the applied in-plane electric field has to 
compete with the spin-orbit coupling induced zero-field 
splitting Ago- Since we expect to find weak coupling 
to electric field and small coupling constant d, it would 
require very strong electric field to achieve the regime 
d|E| ~ Ago- The in-plane components of magnetization 
M x , My, on the other hand show clear dependence on 
electric fields, Fig. [TSJ At low magnetic fields the in- 
plane component of polarization appears and grows with 
the strength of in-plane electric fields. However, the elec- 
tric field dependence becomes less pronounced when an 
additional magnetic field is applied normal to the trian- 
gle's plane. 

In the dependence of the magnetization on electric 
fields, and for the case of an in-plane magnetic field, we 
find the same two main regimes as in the study of the 
polarization: £ k B T (i) and £ <C fc^T (ii). In the first 
case (i) we obtain 



Mi 



ng l n B B l f 1 £Al Q 
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while for the second case (ii) we get 
ngj_fi B B 
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The magnetization shows a strong dependence on the 
electric field E, especially for £ ^> B where this is linear 



in i?-ficld. For low electric fields, however, the magne- 
tization shows only a weak dependence on the electric 
field, both at low and high temperatures. 

For the magnetization (along z) in the presence of a 
perpendicular (also along z) magnetic field we obtain 



M* = ^^tanh(/3B), 



(145) 



which is manifestly independent of the spin-electric cou- 
pling constant d. 



2. Susceptibilities 

The effects of spin-electric coupling on the polariza- 
tion of a molecule show up in the electric susceptibil- 
ity and the spin-electric susceptibility. In Fig. [TB] and 
Fig. [T71 we plot the xx and xy component of the elec- 
tric susceptibility tensor as a function of electric field for 
various strengths and orientations of an additional mag- 
netic field. Both susceptibilities show a high peak in the 
region of weak electric fields that is suppressed by in- 
plane magnetic fields. The peaks are pronounced at low 
temperatures, and vanish as the temperature exceeds the 
splitting of the two lowest-energy levels, k^T ^> d|E|ii. 

In the case of in-plane magnetic field, and weak cou- 
pling to the electric field <i|E| <C Ago, we can calculate 
the electric XEiEj and spin-electric XEiHj susceptibilities 
in the two limiting cases (i) and (ii) defined above, with 
i = x,y. For the electric susceptibility we obtain: 



nd 2 B(£ 2 b %3 - £j£j) 
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XEi Ej 



i so u y 



4A| 



for £ < B 



for £ » B 



(146) 



25 




'0 0.02 0.04 0.06 0.08 0.1 0.02 0.04 0.06 0.08 0.1 

d V A so d V A so 



FIG. 16: Electric susceptibility (xx component), Eq. (I137|l . as 
a function of the electric field in x direction. The three lines 
correspond to various values of the external magnetic field 
in the x direction. The plot is for the temperature feT = 
0.001 Ago- In the inset, the same quantity is plotted at a 
higher temperature, k B T = O.lAso. 



FIG. 18: Linear magnetoelectric tensor (xx component) in 
Eq. (|139p as a function of the electric field in x direction. 
The three lines correspond to various values of the external 
magnetic field in the x direction. The plot is for the temper- 
ature k B T = 0.001 Asci- 




tic. 17: Electric susceptibility (xy component), Eq. (|137|) . as 
a function of the electric field in x direction. The three lines 
correspond to various values of the external magnetic field 
in the x direction. The plot is for the temperature ksT = 
O.OOlAso- In the inset the same quantity is plotted at a 
higher temperature feT = O.lAso- 



in the first case (i), and 



Xe,e 3 



ndA 2 so 
4A| 
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(3A t 
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in the second case (ii). We see that for low E'-fields, 
the electric susceptibility XEiEj depends strongly on the 
applied electric field, and even vanishes if the field is ap- 
plied, say, along x or y directions. For large i?-fields 
instead, the electric susceptibility becomes independent 
of the electric field itself and, for low magnetic fields 
(i.e., for B <C Aso) this reduces to a constant value 
XEiEj = Sijnd 2 /4. At finite (large) temperatures the 
electric susceptibility is still independent of the electric 
field, but it is enhanced by thermal effects ~ 1/T. 

For the electric susceptibilities XEiE in perpendicular 



magnetic field, we obtain 
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where we assumed Aso 3> ^bT, as in the previous Sec- 
tion. As expected, there is no dependence of XEiE 6 on 
the B-field, and for vanishing electric field the electric 
susceptibility reduces to a constant XEiEj = nd 2 /AA^o- 
The quantity of most interest in the present spin sys- 
tem is the spin-electric susceptibility XEiBj , i-e. the mag- 
netic response (electric response) in electric fields (mag- 
netic fields). The nonzero spin-electric susceptibility al- 
lows for the electric control of magnetization and mag- 
netic control of polarization in the crystals of triangular 
MNs, even in the case when the coupling between the 
molecules is negligible. In addition, XEiBj is nonzero 
only in the presence of spin-electric coupling, i.e. when 
d + 0. 

The spin-electric susceptibility shows a characteristic 
peak in weak electric fields which vanishes in an exter- 
nal magnetic field, see Figs. [18] and [19] The peak in 
the diagonal ^-component, Xe x e w , moves towards the 
higher electric fields and broadens as the magnetic field 
B x increases. The peak in the off-diagonal component 
Xe x b z , on the other hand, shifts towards the lower elec- 
tric fields, and narrows as the in-plane magnetic field 
increases. Both peaks disappear at high temperatures, 
/c B T> A so . 

For in-plane magnetic fields and weak spin-electric cou- 
pling the spin-electric susceptibility XEiBj is 
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for the low temperature case (i) , while for the second case 
(ii) we obtain 
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FIG. 19: Linear magnetoelectric tensor (xz component) in 
Eq. (|139[) as a function of the electric field in x direction. 
The three lines correspond to various values of the external 
electric field in the x direction. The plot is for the temperature 
k B T = O.OOlAso. 



FIG. 21: Magnetic susceptibility (zz component) in Eq. 
(|138[) as a function of the magnetic field in x direction. The 
three lines correspond to various values of the external elec- 
tric field in the x direction. The plot is for the temperature 
feT = O.OOlAso- The inset represents the same quantity at 
a higher temnerature knT = 0.1 A«n. 
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FIG. 20: Linear magnetoelectric tensor (xx component) in 
Eq. (|139p as a function of the magnetic field in x direction. 
The three lines correspond to various values of the external 
electric field in the x direction. The plot is for the temperature 
k B T = O.OOlAso- 



FIG. 22: Magnetic susceptibility (xx component), Eq. (|138p 
as a function of the electric field in x direction. The three 
lines correspond to various values of the additional magnetic 
field in the z direction. The plot is for the temperature feT = 
O.OOlAso- 



By inspecting the above expression, we can infer that fol- 
low temperatures and low .E-fields the spin-electric sus- 
ceptibility shows no dependence on the absolute value of 
the electric field E and only a weak dependence on the 
applied magnetic field B. Moreover, when both fields are 
applied along one special direction, say, along x, and as- 
suming also B <C Aso, the spin-electric susceptibility be- 
comes xe x B m = ndgi/iB /4Aso> i-e. it reaches a constant 
value. The finite temperature expression shows that the 
spin-electric response is reduced, as opposed to the elec- 
tric response where temperature increases the response. 
Thus, for strong spin-electric response one should probe 
the spin system at low temperatures (fc^T <C Ago)- 

The diagonal out-of-plane component of the magnetic 
susceptibility, xs^s^in the presence of an external mag- 
netic field in the x direction decays strongly in the ap- 
plied electric field along the x direction, Fig. [5T] In 
electric fields, the xs x ,b x component shows a peak that 
is reduced by the application of the magnetic field in x 



direction, Fig. 

We can derive the magnetic susceptibilities in the two 
regimes. In the first case (i) we obtain (assuming now 
only linear effects in E- field): 
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with i,j = x,y, while 
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for S z Ago <C BE. At low temperatures the in-plane mag- 
netic susceptibility shows a linear dependence on the ap- 
plied electric field E, thus allowing for a simple estimate 
of the electric dipole parameter d from magnetic mea- 
surements. Note that for strong electric fields (8 ^> B), 
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the magnetic susceptibility can vanish, since the mag- 
netization does not depend on the magnetic field any- 
more. However, such a regime would not help to identify 
the electric dipole coupling strength d from susceptibil- 
ity measurements directly. The perpendicular magnetic 
susceptibility shows a strong electric field dependence 
XBzBz ~ £~ x and can be used as an efficient probe for 
extracting the electric dipole parameter d. In the second 
case (ii) we obtain 
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FIG. 23: Heat capacity, Eq. (|157|) . as a function of tempera- 
ture in various electric fields. 



when i, j = x, y, and 
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We consider the cases of perpendicular B-field and in- 
(154) plane B-field in the limit Aso 3> k B T. In the first case, 
i.e. for B II z we obtain 



The magnetic response increases with temperature. Also, 
in this limit the dependence of the magnetic susceptibility 
on the applied electric field is rather weak (xBiBj{£) ~ 
E 2 ), thus this regime is also not suitable for observing 
spin-electric effects. 

For the magnetic susceptibility in the perpendicular 
magnetic field we find 



XB Z B Z 
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while for the in-plane magnetic susceptibility xb i(s) b x(j) 
we obtain 
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in the limit B,k B T <C Aso- We mention that for B 
perpendicular to the molecular plane there is no electric 
field E (magnetic field B) dependence of the magnetiza- 
tion Mi (electric polarization Pi). Thus, in order to see 
spin-electric effects one needs to apply magnetic fields 
which have non-zero in-plane components. 



D. Heat capacity 

Next we investigate the dependence of the heat ca- 
pacity on the applied electric and magnetic fields in dif- 
ferent regimes. The heat capacity is defined as C — 
-d/dT(dhx(Z)/d0), so that we obtain 
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The heat capacity C shows a quadratic dependence on 
the applied electric field for the entire range of E-field 
strengths. On the other hand, the magnetic field depen- 
dence of C is non-monotonic, and shows a maximum for 
some finite B-field strength B max ~ fc^T. In the second 
situation, i.e. for Biz we get 



C ~ Nk B f3 2 
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£ < k B T. 



As in the previous case, the dependence of the heat- 
capacity C is linear in E-field for low E-fields. How- 
ever, for large E-fields the dependence is non-monotonic 
and thus shows a maximum for some finite electric field 
strength £ max ~ k B T. Note that in this case also the 
dependence on the magnetic field is non-monotonic, and 
thus we obtain a second maximum for B max ~ k B T. 
We can conclude from the above expressions that the 
strongest dependence of the heat capacity on the electric 
field is when the magnetic field is applied in-plane, and 
then it is mostly quadratic. 

For the derivation of all the thermodynamic quanti- 
ties presented in the previous sections, we have restricted 
ourselves to the contributions arising from only the low- 
est four states, even though the spin system spans eight 
states in total. This description is valid if the splitting 
between the energies of S — 1/2 and S — 3/2 states 
is much larger than the temperature kBT. This split- 
ting varies strongly with the applied magnetic field, for 
B = 3J/4 one of the S = 3/2 states (M = -3/2) cross- 
ing the M = 1/2 of the 5 = 1/2 states and, even more, 
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FIG. 24: Heat capacity, Eq. (|157[) . at low temperature as a 
function of external electric field. 



for B > 3 J/2 the M = —3/2 becomes the spin system 
ground states. Thus, for large magnetic fields our effec- 
tive description in terms of only the 5 = 1/2 states breaks 
down and one have to reconsider the previous quantities 
in this limit. 



V. CONCLUSIONS 

Electric fields can be applied at very short spatial and 
temporal scales which makes them preferable for quan- 
tum information processing applications over the more 
standard magnetic fields . Nanoscale magnets, while dis- 
playing rich quantum dynamics, have not yet been shown 
to respond to electric fields in experiments. We have 
investigated theoretically the possibility of spin-electric 
coupling in nanoscale magnets using symmetry analy- 
sis, and found that the spin-electric coupling is possible 
in antiferromagnetic ground-state manifolds of spin- 1/2 
and spin-3/2 triangles, as well as in spin-1/2 pentagon. 
The spin-electric coupling in the triangle can exist even 
in the absence of spin-orbit coupling, while the coupling 
in the pentagon requires the spin-orbit interaction in the 
molecule. We have characterized the form of the spin- 
electric coupling in all of these molecules and presented 
the selection rules for the transitions between the spin 
states induced by electric fields. 

While the symmetry can predict the presence or ab- 
sence of the spin-electric coupling, it can not predict the 
size of the corresponding coupling constant. In order to 
find a molecule suitable for electric manipulation, it is 
necessary to have an estimate of the spin-electric cou- 
pling strength. For this purpose, we have described the 
nanoscale magnets in terms of the Hubbard model, and 
related the coupling constants of the symmetry-based 
models to the hopping and on-site energy parameters of 
the Hubbard model. We have found that the modifica- 
tion of the Hubbard model parameters due to the electric 
field produces a spin-electric coupling of the same form as 
predicted by the symmetry analysis. However, within the 
Hubbard model, the coupling constants have a clear and 



intuitive meaning in terms of the hopping and on-site en- 
ergies of the localized electrons. We have also studied the 
superexchange interaction of the spins on the magnetic 
centers through the bridge. If we assume that the inter- 
action of the localized spins is a property of the bridge 
alone, the spin-electric coupling can be calculated by ab- 
initio analysis of the bridge alone, and not of the entire 
molecule. 

Finally, we analyzed the role of spin-electric coupling 
in standard experimental setups typically used for the 
characterization of nanoscale magnets. We find that the 
spin-electric coupling can be detected in the ESR and 
NMR spectra that probe the local spins. Also, ther- 
modynamic quantities, like the polarization, magnetiza- 
tion, linear magnetoelectric effect, and the specific heat 
show signatures of spin-electric coupling in the triangular 
molecules. Thus, our results set a path toward finding 
suitable molecules that exhibit spin-electric effects and 
how they can be identified experimentally. 

In this work, we have focused on the spin rings with an 
odd number of magnetic centers (odd spin rings), whose 
low-energy spectrum is dominated by frustration effects. 
The odd spin rings, due to frustration, possess a four-fold 
degenerate ground state multiplet, which can be split 
by electric fields. As opposed to the odd spin rings, 
the ground states of even-spin rings is usually a non- 
degenerate S = state, separated from the higher energy 
states by a gap of the order of exchange coupling J. Cou- 
pling of the electric field to these states can thus proceed 
only via excited states, and the coupling strength is re- 
duced by d\Ei\/ J. Similarly, in lower-symmetry odd-spin 
rings, the ground state multiplet consists of an S = 1/2 
Kramers doublet, which can not be split by electric fields, 
i.e. there is no spin-electric effect in zero magnetic field. 
Therefore, the odd spin rings seem to be the most suit- 
able candidates for observing the spin-electric coupling 
and using it to control the spins. 
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sions. We acknowledge financial support from the Swiss 
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Appendix A: Spin states in terms of the c r operators 



In this appendix we show the expressions for the three- 
electron symmetry adapted states iV'r' 7 ) m EcjS- (|35|) and 
(|38p in terms of the symmetry adapted creation operators 
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cj, a . Making use of Eq. (|47l) we obtain 
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where a stands above for sign(a). 

Appendix B: Hso, H°_^, and H^_ d matrix elements 

For the SOI matrix elements we obtain 



while the remaining terms are equal to zero. For the 
electric dipole matrix elements we obtain 
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